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Abstract 

Qh, The Monte Carlo program KoralW version 1.42 is presented. It generates all 

Q^l four-fermion final states with multibranch dedicated Monte Carlo pre-samplers 

^ ' and complete, massive, Born matrix elements. The presamplers cover the entire 

phase space. Multiphoton bremsstrahlung is implemented in the ISR approximation 
within the YFS formulation with the 0{a^) leading-log matrix element. The anoma- 
/\ ' lous WWV couplings are implemented in CC03 approximation. The standard de- 

» . cay libraries (JETSET, PHOTOS, TAUOLA) are interfaced. The semi-analytical 

CC03-type code KorWan for diff'erential and total cross-sections is included. 
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NEW VERSION SUMMARY 

Title of the program: KoralW, version 1.42. 

Reference to original program: Comput. Phys. Commun. 94 (1996) 215. 
Authors of original program: M. Skrzypek, S. Jadach, W. Placzek and Z. Wcis 
Computer: any computer with the FORTRAN 77 compiler and the UNIX operating 
system 

Operating system: UNIX (program tested under AIX 4.x, HP-UX 10.x), Linux 
Programming language used: FORTRAN 77 
High-speed storage required: < 10 MB 

No. of cards in combined program and test deck: about 30,000 plus 11,261+5,970+2,264 of 
physics generators libraries (JETSET+TAUOLA+PHOTOS) plus 369,331 of an external 
matrix element library (GRACE). 

Keywords: Radiative corrections, initial-state radiation (ISR), leading-logarithms (LL) 
approximation, heavy boson W, 4-fermion processes, Monte Carlo (MC) simulation/genera- 
tion, quantum electrodynamics (QED), quantum chromodynamics (QCD), Yennie-Frau- 
tschi-Suura (YFS) exponentiation. Standard Model (SM), LEP2. 

Nature of the physical problem: The W^-pair production and decay is and will be used 
as an important data point for precise tests of the standard electroweak theory at LEP2 
and higher energies. The effects due to background processes, QED bremsstrahlung and 
apparatus efficiency have to be subtracted from the data. The program deals with all 
e~^e~ processes leading to 4-fermion final states accompanied with multiphoton initial- 
state radiation. It also includes the effects of the Coulomb correction, 'naive' QCD, 
anomalous couplings, quarks hadronization, r decays, and photon radiation in leptonic 
decays. 

Method of solution: The Monte Carlo methods are used to simulate all 4-fermion final- 
state processes in the e~^e~ collisions in the presence of multiphoton initial-state radiation. 
The latter is described in the framework of the YFS exclusive exponentiation. The Im- 
pair production is included in a 'natural' way as a subset of the Feynman diagrams for 
the above processes, but it can also be generated exclusively by switching to the so-called 
CC03 process. The Monte Carlo generation is done on an event-by-event basis, with 
constant or variable weights, where an event is represented by flavours and four-momenta 
of all respective particles - supplemented with a collection of weights, if the variable 
weight option is chosen. After the event generation is completed the program provides 
the cross sections together with their statistical errors for all the processes involved. Any 
experimental cuts and apparatus efficiencies may be introduced easily by rejecting some 
of the generated events. 

Restrictions on the complexity of the problem: Only processes with 4-fermion final states 
are considered. QED radiative corrections are implemented in terms of multiphoton ISR 
in the YFS Monte Carlo framework with the 0{a^) LL-type matrix element. For the 
CC03 subset of diagrams the Coulomb correction for the intermediate WW states is 
also included. The final-state QED radiation is generated for charged leptons with the 
help of the program PHOTOS in the LL approximation (up to two photons). QCD 



effects are included in the so-called "naive QCD" approximation. A part of electroweak 
corrections is incorporated in the "improved Born approximation" (through appropriate 
renormalization scheme). Anomalous triple gauge boson couplings are implemented only 
in the CC03 subset of diagrams (i.e. the W^-pair production). Quadruple gauge boson 
couplings are not implemented. The r-decays and quark hadronization are performed, 
respectively, with the help of the dedicated packages TAUOLA and JETSET. 
Typical running tinae: On IBM PowerPC M43P240 (266 MHz, 65 CERN units) in- 
stallation one needs: (a) 2.5 sec per 1000 constant-weight events and 0.6 sec per 1000 
variable-weight events for CC03 matrix element and (b) 12500 sec per 1000 constant- 
weight events and 6 sec per 1000 variable-weight events for a complete 4-fermion matrix 
element (GRACE). These results are for a default/recommended setting of input param- 
eters but with all decay libraries switched OFF. 



1 Introduction 

As LEP2 approaches the 200 GeV centre-of-inass energy regime, the need for rehable 
calculations of the so-called signal CC03 and general background processes in e~^e~ — >■ 
W^W~ —>■ fif^fsJA becomes more and more immediate. Moreover, with the added need 
for the accommodation of arbitrary detector cuts in these calculations, the only practi- 
cal solution is the Monte Carlo event generator realization of the calculation, wherein 
these cuts may be imposed on an event-by-event basis. In the present paper, we provide, 
in version 1.42 of the program KoralW [^ such an event generator in which all rele- 
vant A-fermion processes, both charged-current (CC) and neutral-current (NC) ones, in 
e~^e~ — i> /1/2/3/4, in the LEP2 200 GeV regime are realized in the presence of the YFS @ 
exponentiated initial-state multiple photon radiation (ISR). 

More precisely, in the WW sector this version of the program has as its ultimate 
objective the sub-per-cent precision regime of 0.5% as called for in the LEP2 Workshop 
Report of Ref. [Q. Thus, when compared with version 1.02 [|l|, the inclusion of the 4- 
fermion background processes for the e"'"e~ — > W^W~ — > /1/2/3/4 signal is an essential 
improvement. This is described in detail in Section 3. In addition, new physics anoma- 
lous couplings are now featured for the CC03 class of graphs - these were absent from 
version 1.02. The ISR is now calculated through the 0{a^) LL in the YFS exponenti- 
ated framework; in Ref. [^, it was only calculated through the O(a^) LL in the same 
framework. The Coulomb correction is now implemented in a way that is more reliable 
near the threshold, following the authors of Ref. Qj. The leading QCD correction, as 
well as the CKM matrix, is now featured following Ref. 0; in version 1.02, these correc- 
tions are absent. Two other effects featured in the current version, which are missing in 
version 1.02, are the colour reconnection effect, which we model after Ref. P], and the 
Bose-Einstein effect, which we take after the authors of Ref. 0. In addition, three choices 
of renormalization scheme are available. Considerable technical precision checks of the 
program have been made, using semi-analytical results as described below, so that we 
have established the technical precision of the program at the level of 0.2%. The primary 
missing ingredient for the final step to a total precision tag of 0.5% is the implementation 
of the exact 0{a) electroweak corrections for the CC03 class of graphs in the YFS expo- 
nentiated framework, with the consequent intermediate-state n{'-f) soft radiative effects, 
as we have already published using the program YFSWW3-l.il in Ref. P|. These effects 
will be incorporated in a later version of KoralW |^ . As a result of the studies done in 
Refs. [|, |1^, and for the preparation the current paper, we can already safely set, in 
the current version 1.42 of KoralW, the precision tag for the cross section normalization, 
for the process e^e^ -^ W~^W~ — > /1/2/3/4, to 2% in the LEP2 energy regime, where 
it is understood that a cut of some type is used to define the two W^s. The same 2% 
precision tag should remain valid also for the Z-Z physics, even though a formal study 
is missing. This, however, may not be true in general 4-fermion final-state processes. 
Here, the precision may depend very much on a particular final state under consideration, 
a choice of physical observables, experimental cuts, events selection criteria, etc. Thus, it 
requires a more dedicated study. The most problematic are the processes with electron(s) 



and/or positron(s) in the final state, particularly when regions of small scattering angles 
are allowed (e.g. selection criteria allow some of the particles to be lost in the beam 
pipe) and/or events with high-p^ photons are accepted. Various problems that can be en- 
countered there are discussed in some detail in Subsections p.2.4| and p.3.2| . Nevertheless, 



at the Born level, our program can be used in all regions of the phase space, including 
the most singular ones. Also, an important technical step has been achieved, which is 
indispensable for a reliable solution of bremsstrahlung implementation in the general case. 

Let us also point out that contributions of the A-fermion processes from some regions 
of the phase space are implicitly included in 2-fermion generators as a part of radiative 
corrections (analytically cancelled out with the leading (9(q;^) virtual corrections). One 
needs to be cautious about this point so as to avoid double counting. 

The outline of the paper is as follows. We present in detail the overall structure of the 
program in Section 2. In Section 3, we fully describe the 4-fermion phase-space genera- 
tion and the corresponding matrix elements used in the calculations of the program, the 
generation of the energy distribution due to ISR in the program, and the semi-analytical 
CC03 distributions contained in the program. Section 4 describes the practical use of the 
program, so that it should help the user to take advantage of the program's capabilities 
in an efficient and easy manner. Appendices contain useful technical information on the 
construction and use of the program, its matrix elements, its input/output, etc. 

2 Structure of the Program 

In this section we provide the reader with a brief guide of the KoralW program. We 
will describe its main routines, libraries and interfaces. We want to note here that the 
program, in its distribution version, is prepared for a UNIX-type operating system that 
supports directories and make command. Namely, we have divided the source code into 
a number of subdirectories, in order to make the structure more transparent and easier 
to handle. Also the structure of Makef ile-s is provided, for easier compilation, as well 
as some other auxiliary functions (e.g. clean-up). This structure, which is in fact not 
very complicated, can be in principle avoided and code can be transformed into a single 
FORTRAN file. Some care must be taken while handling the ampli4f .grc.all routines 
where certain include files with the same name, have different contents in different sub- 
directories. Also, note that the following three directories: ampli4f, ampli4f .grc.all 
and ampli4f .grc. all-old are mutually exclusive, i.e. they contain three versions of the 
same routines and only one at a time can be used by the program. 

2.1 ainpli4f — Template of the Library 

Dummy directory for the external matrix element library, to be filled in by the user. 
The front-end routines ampini and amp4f are expected by KoralW to be put here by the 
user. We expect the user to replace this directory with his/her own code, e.g. with their 
favourite parametrization of anomalous couplings. 



2.2 ampli4f .grc.all — All Four-Fermion Library 

An implemented library of external matrix elements for all four-fermion final states has 
been generated by the GRACE v. 2 package [0]. This dedicated code has been provided 



for KoralW with the courtesy of the Minami-Tateya Group of KEK. 

The directory amp4f contains routines for the actual computation of matrix elements. 
The auxiliary library is located in the directory channel. The routines specific to KoralW, 
which transmit and reset parameters and dip-switches of GRACE routines are located in 
the directory grc4f _init. Finally, the front-end routines ampini and amp4f are in the 
amp4f directory. 

2.3 ampli4f . grc . all-old — Old CC-all Library 

An old version^ of the library of external matrix elements for all lypF-type final states. 



generated by the GRACE v. 1 package [|T3| . This library contains all the scripts capable of 
building the FORTRAN code from scratch with the help of the symbolic package GRACE. 
The scripts are located in the subdirectory grace4f _init. All of the FORTRAN files and 
directories created by the scripts can be identified by the word expression ".auto." in 
their names. Similarly, the template files used by the scripts are marked with the word 
expression ".template.". 

2.4 demo — Demonstration Programs 

This subdirectory contains two demonstration subprograms: KWdemo and KWdemo2. Gen- 
erally the user is supposed to provide his/her own main program. Nevertheless the two 
subprograms quoted here are simple examples of a main program. 

KWdemo is the first of them. It has a double role as a useful template and as a first 
cross-check that the MC generator KoralW runs correctly on a given installation. The 
essential part of this program is a loop in which a series of KoralW events is generated. 
It also reads the input from a disk file, but no histogramming is performed and most 
of output is from the generator itself. At the end of the program, a MC integrated 
cross section of KoralW is compared with a semi-analytical result from KorWan. The 
program is compiled/linked and executed, with the help of the makefile, for two data sets, 
as follows: maike KWdemoCCOS and make KWdemoGRCall. The two data sets for the two 
separate runs are: 

• demo . 14x/190gev/KWdemo . input .CC03 for which the CC03 "signal" process is sim- 
ulated (internal matrix element), the ISR is on, with hadronization, tau decays, no 
CKM mixing, weighted events. 

• KWdemo . input . GRCall for all four-fermion matrix elements form GRACE (external 
ME), hadronization, tau decays, weighted events. 



See Sect. |3.2.3 for explanation why we have decided to keep this "old" package. 



KWdemo2 is the second example. Apart from the loop in which series of KoralW events 
is produced the example of histogramming with proper normalization is also included. It 
is done in the subprogram ROBOL. The program is compiled/linked and executed, with 
the help of the makefile, for two data sets, as follows: make KWdemo2HADR and make 
KWdemo2SEMI. The two data sets provided for the two separate runs are: 

• KWdemo2 . input . HADR CC03 "signal" process is simulated, (internal matrix element), 
the ISR is on, all four-quark channels only - selected with Umask matrix (see Ap- 
pendix B for a definition), WtMain = 1 events. 

• KWdemo2 . input . SEMI CC03 "signal" process is simulated (internal matrix element), 
the ISR is on, anomalous triple gauge boson coupling constants are activated, see 
data file, semi-leptonic channel only - selected with the Umask matrix. Variable- 
weight events are generated and analysed. 

2.5 glib — Histogramming 

A handy FORTRAN histogramming package, GlibK ||T^, is provided in this directory. 
It is used by KoralW both for hard-coded internal bookkeeping and for some optional 
"external" tests. The package is similar in its usage to the classic HBOOK package of 
CERNLIB. 

2.6 interfaces — Interfaces to the Libraries 

Interfaces to external libraries are collected in this directory. 

1. tohep fills in the /HEPEVT/ common block, decays r leptons (with TAUOLA) and 
generates bremsstrahlung in W decays as well as secondary decays (with PHOTOS). 

2. tohad does the hadronization (with the help of JETSET) along with necessary 
colour (re) connection. 

3. The interface to the external matrix elements is provided by the routines ampinw and 
ampext. The first performs the necessary initializations. It calls the user-supplied 
initialization routine ampini. Some additional parameters (apart from xpar) can 
be reached by the user with the help of masow and kwparm2 routines. The ampext 
transmits to KoralW the value of the external matrix element that it calculates 
by a call to the user-supplied routine amp4f . It also allows, with the help of the 
dip-switch key_cms_eff , choosing the four- vectors supplied for this calculation to 
be in the LAB or effective CMS (default) frames. 

2.7 kwlund— JETSET and PHOTOS libraries 

JETSET V. 7.4 [J15l and PHOTOS [|l^] are located here. Note that the common block 



/HEPEVT/ is expected to contain single-precision (REAL*4) variables. 
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2.8 model — Matrix Elements 

The "model" weights are calculated here. 

1. The CC03 matrix element is calculated by wwborn. It takes as an input massive 
four- vectors of final fermions (assumed to be in their rest frame with beams along the 
z-axis, e~ in the +z direction). Next, it "reduces" these four- vectors to the mass- 
less one^ in three ways. The "sophisticated" method, which reconstructs angles 
(invkin) and rebuilds four- vectors with the same angles and zero masses (kineww). 
The next method simply rescales the three-momenta, thus breaking the overall mo- 
mentum conservation. The third provided option is "no reduction at all". These 
massless four-vectors are transmitted to wwborn_massless that does the actual cal- 
culation with the help of the wwprod and wdecay routines. If anomalous couplings 
are requested the wwamgc routine is called instead of wwprod. 

2. The ISR photonic corrections (the so-called beta-functions) up to the third order are 
provided by the routine betar, which in turn calls the d_isr* routines to construct 
the actual real and virtual bremsstrahlung contributions. 

3. The Coulomb correction is located in culmc routine. It is used by Monte Carlo 
routines. Semi-analytical routines have an identical correction implemented in the 
culsan routine in semian directory. 

2.9 semian — Semi-analytical Routines 

This directory contains a package for semi-analytical calculations. 

1. The total cross section a, with or without ISR, is provided by the korwan routine. 
For the differential da/d\ogv it uses the yfspho function. For the total a with 
ISR, korwan integrates yfspho and, in the case of no radiation, a is provided by 
the function xsmuta; this in turn integrates the dlmuta, which integrates d2muta, 
the actual two-dimensional differential cross-section. The function yfspho provides 
various kinds of one- dimensional photonic distributions da/dv (with soft residual 
subtracted) or da/dlogv. It uses xsmuta as well. 

2. The one-dimensional distribution of the single- W^ invariant mass is provided by the 
slwan function. 

3. The two-dimensional distribution of the double- TV invariant mass is provided by the 
sls2wan function. 

4. The average mass, (l/o") j dvdsids2{^/si + y/s2 — 2Mw)da /{dvdsids2), is calculated 
by the mavrg routine with the help of korwan with negative s- variable input param- 
eter. 



Note that for the CC-all and NC-aU modes the fully massive matrix element is used, see Sect. p. 2.2 



5. The average mass loss, {^/s|2){l/a) J dv v da/dv, is calculated by the mloss routine 
with the help of korwan with negative keypho input parameter. 

2.10 tauola — TAUOLA library 

It is the directory with the TAUOLA library |T^ for simulation of r-decays. Note that, 
as usual, in the distribution version of TAUOLA the parameters in the r-decay modes are 
not adjusted appropriately. We recommend that user replaces this version of TAUOLA 
by the version of his/her own collaboration]^. 

2.11 korww — "Central Processing Unit" 

This directory contains the actual Monte Carlo event generator. Subprograms used di- 
rectly by the user are the following: 

• KW_ReaDataX - the subprogram used to read, from the disk file, the default 
input data of KoralW and subsequently the data of the user into the array xpar 
at the very beginning of the use of KoralW. (This subprogram did not exist in the 
version 1.41.) 

• KWJnitialize - the subprogram that does all initializations of internal variables 
- it calls several initializers of the main internal modules of the generator, such 
as karlud(-l, . . .), and of TAUOLA and the external matrix element. It prints 
out directly or indirectly all the input parameters. (This subprogram replaces the 
koralw(-l,xpar,npar) entry in version 1.41.) 

• KW_Make - the most important subprogram of KoralW. It generates single MC 
events. Functionally it is the high-level management subprogram in the event gen- 
eration. It invokes: 

1. karlud, which provides the four-momenta of the outgoing fermions and the 
ISR photons along with the value of the crude distribution, 

2. KW_model_4f , which computes the Born matrix elements either CC03 or CC/NC- 
all, and all additional effects requested by the user, e.g. anomalous couplings 
and/or the Coulomb correction, 

3. betar, which calculates the QED ISR corrections up to third order, 

4. tohep, tohad, the programs from TAUOLA, PHOTOS and JETSET libraries, 
which perform, respectively, r-decays, generation of bremsstrahlung in W and 
r decays, and hadronization. 

KW_Make also decides about the rejection of an event in the case of a constant- weight 
event, or of a semi-constant-weight event. 



In case it is impossible, please contact the authors. 



KW_Finalize does all final bookkeeping, including the calculation of the integrated 
(total) cross section. It prints a summary output on all series of generated MC 
events. (This subprogram replaces the koralwd ,xpar ,npar) entry in version 1.41.) 



The other important internal units of the generator are: 



karlud The routine karlud manages the actual phase-space point generation. Upon 
initialization in mode=— 1, the event generation sequence in mode=0 is the following: 

1. yf sgen generates the s' variable and the ISR photons in the LAB frame. 

2. decay generates the decay channel based on pretabulated cross-sections (and 
stores its ID number in the variable label). 

3. make_phsp_point_* generates four-vectors of fermions depending on the type 
of the final states drawn. It takes place in the effective CMS frame, i.e. the rest 
frame of the outgoing fermions. There are two different generators provided in 
this place. 

4. from_cms_ef f transforms the four-momenta of the event from the effective 
CMS to the LAB frame. 

5. selecto imposes cuts on the four-momenta. The selecto routine contains 
the built-in cut-offs, whereas the user_selecto routine is provided for the 
user-defined cuts - it is located in the file demo. 14x/user_selecto.f and by 
default no cuts are applied there. It is a matter of efficiency that as many of the 
cuts as possible should be imposed in this low-level routine to avoid further, 
time-consuming, steps of event generation for events that are to be rejected. 

6. get_phsp_weight_* calculates the crude distribution (the crude weight) for the 
accepted event. 

7. karlud calculates necessary overall normalization factors. 

inode= 1 and 2 provide, as usual, some post-generation information. 
In the following * is a wildcard for a part of a routine name. 

make_phsp_point_* calls the *_brancher routines that do the random choice of 
kinematical branch to be generated according to the preset probabilities. It should 
be noted that these probabilities are dummy parameters and can be changed by the 
user, e.g. to cross-check the program or to fine-tune its efficiency. It is also possible 
to add some more branches to the generator. Next, make_phsp_point_* calls the 
*_spacegen routines in order to build the actual four-momenta. 

get_phsp_weight_* performs the summation over various branches of the multi- 
branch four-fermion crude distribution in order to construct the total normalization 
of the event. It also uses the *_spacegen routines. 



• *_spacegen. The actual hard-core of the generation and normahzation of the 
four-fermion phase-space point. It is based on a comphcated series of kinematical- 
variables generation. It is performed in various frames, in various orders and with 
different types of singular behaviour, depending on the final state and, subsequently, 
on the branch chosen at a given time by the program on previous stages of event 
generation. Details of the algorithm will be presented elsewhere . 

2.12 B.E. directory for the Bose-Einstein effect 

The package of programs in subdirectory B.E. implements the Bose-Einstein (BE) effect 
in the hadronization using, the "weight method" described in Ref. [0. It is not the 
integral part of the MC generator but rather a stand-alone application of KoralW. It has 
a double role: (1) to illustrate how to implement the Bose-Einstein effect according to the 
method of Ref. [0, (2) to provide an example (template) of the use of KoralW in a C++ 
environment. The entire exercise on BE is implemented in C-I--I-. KoralW provides "raw" 
events, which are translated to C++ structures, and the whole process of constructing 
the BE weight, analysing events (i.e. defining jets, eliminating combinatorial background, 
fitting the W mass), histogramming and graphical output is done in C++. There are 
only several lines of Fortran code (in ReaData.f) in the entire subdirectory B.E. The 
histogramming, fitting and graphics are done with the novel ROOT system |T^, also 
entirely in C++. 

The topography of the B.E. directory is the following: 

B.E./src - sources and compilation/link objects, 

B.E. /run - data files and outputs for /from runs, 

B.E. /fig - analysing results from . ./run, all plotting/fitting, etc., 

B.E. /bak - attics. 
See also README file in B.E. directory. 

The source files and executables are in B.E./src directory. Let us very briefly char- 
acterize their functionality: 

• rmain . C is the main program. It runs the main loop over events. 

• Semaph class is used by rmain . C to manage main loop over events using information 
from "semaphore" disk filesQ. 

• ROBOL class procedures are called in rmain . C. It manages the generation of events, 
the calculation of the BE weight, and the analysis of MC events, all by dedicated 
classes. 

• KoralwMaker class provides for the C++/F77 interface with the KoralW event gen- 
erator. It transfers input data from the main program, which is in C++, to the 
KoralW generator, which is in F77. It also picks up the output data (events) from 
KoralW and makes them accessible to ROBOL and other C++ classes. 



*The Semaph class is a translation from the analogous F77 subprogram in BHLUMI p9| ]. 
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• KorEvent is a class for the single MC event from KoralW, with some additional 
data fields/members for the analysis. The ROBOL procedures handle the current MC 
event as a KorEvent object. 

• PartLund is a class for a single particle in a format very close to the PDG/LUND 
convention - a single event of the type KorEvent is simply a list of objects of the 
type PartLund. 

• VLorentz is a class for describing the single Lorentz four-momentum. It provides 
some basic functionality such as linear algebra and dot-product. It is used to con- 
struct the type PartLund in other places of the package. 

• BEwtMaker is a class in which the current event defined by the KoralwMaker object 
gets assigned the BE weight. 

• JetAnalyzer is a class in which the current event undergoes standard analysis - 
four jets are defined, the combinatorial background is eliminated and the masses of 
jet pairs are recorded in histograms/ntuples. 

Note that F77 externals are referred to only in KoralwMaker and Semaph classes. The only 
extra F77 code is composed of two tiny routines in ReaData. f (they interface OPEN/CLOSE 
functions of F77). 

The input data for the run and the resulting histograms are stored in run subdirectory. 
In particular the run/172GeV.4J is the most important subdirectory. 

The final analysis of the MC results is done in fig subdirectory, using several C++/CINT 
scripts processed by ROOT. In particular all eight figures of Ref. are produced by 
them. For instance the plot with the fit of the two-jet distribution with the Breit- 
Wigner formula is made by the macro view-Fit JetBE.C. It reads histograms from the 
run/200GeV . 4 J/rmain . root file, which can be produced with "make 172GeV . 4 J-start" . 
For more details on how to run the pro duct ion- analysis sequence of the programs, see the 
files README and src/Makef ile. Let us finally advertise that it is possible, in particular, 
to get a menu of graphics programs with the command "make fig" and to produce the 
on-line documentation of all our C++ classes in the html format and view it with the 
html browser with a single command: "make dok-view" . 

In the preparation on the F77/C++ interface we profited a lot from inspecting the 
C++ code of ATLFAST [^ and from discussions with the authors of this program. 
Helpful discussions with the authors of the ROOT system are also acknowledged. 

3 Details of the Program 

3.1 Four-Fermion Phase-Space Generation 

The generation of the four-fermion phase space in KoralW is based on a multibranch type 
of Monte Carlo algorithm, cf. e.g. |[T7|,plH2^- The cross-section is calculated in the usual 
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way, as an average of the ratio of the exact matrix element to the crude one, averaged 
over the crude distribution: 

a = /"rfPhsp|M|2 = /^\ [dp, (1) 

J \ fcR I dp J 

branches 

dp = dPhspfcR= Yl dFhsp' Pi fcR, (2) 

i 

where pi is a probabihty of generating the branch i and dPhsp* is aheady parametrized 
by angles and masses defined in branch-dependent Lorentz frames: 

dFhsp' = ds\dsi JJ d cos 6] dcj)] \]. (3) 

i=l,2,3 

Equation (H) describes a general framework for generating arbitrary phase-space config- 
urations. The process specific information (and the difficulty!) is hidden in the /^^ 
functions. Each branch of /^^ is designed to describe a certain type of singularity that 
is encountered in the Feynman graphs. As an example, the t-channel singularities can 
be thought of as 1/t, as the s-channel ones as resonances, etc. But this is of course 
process-dependent. In the case of four-fermion final states there is nearly a hundred dif- 
ferent possible final states, each of them having up to over a hundred Feynman graphs 
in the matrix element, each graph having a different structure of singularities. Of course 
there is a lot of symmetries and similarities that can be employed in constructing the 
generator branches. Nonetheless, the number of different branches exceeds fifty in the 
case of KoralW. In fact each branch i of Eq. (0) consists of many "sub-branches" , which 
split on subsequent lower levels of generation. The total number of such "elementary" 
branches would exceed a million. Just by looking at this large number, one can realize 
that optimalization must play an essential role in the algorithm. Indeed, there is a large 
number of internal parameters and coefficients, as pi of Eq. (^, that must be fine-tuned 
in order to enhance the desired singularities and damp the others for given final-state 
configurations. 

In the KoralW code we have implemented two different sets of libraries of these f}jp^ 
functions. In addition there is a third possibility of using both libraries simultaneously 
in a stochastic mixture. These packages provide a powerful tool for controlling the com- 
plicated and sometimes numerically unstable integration. Moreover, the coefficients pi of 
Eq. (0) are dummy parameters of the generators and can be changed by the user for the 
consistency tests. 

Details of the actual algorithms will be presented elsewhere 0. 

3.2 External Matrix Element 

In its present version, KoralW includes an interface to the external library calculating the 
correction weight due to a different, external, matrix element. The idea behind this is that 
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the user may occasionally wish to replace the internal matrix element by a different one, 
for instance including special combinations of the anomalous couplings. Thanks to the 
modular structure of KoralW and, in particular, to the factorizability]^ of the approximate 
QED matrix element into the Born matrix element and the QED part, it is straightforward 
to replace the existing Born-level matrix element with any other one, provided that the 
external library is able to calculate the corresponding matrix elements out of the externally 
generated four-momenta. To this end the external program provided by the user must be 
able to calculate the matrix element, completely normalized (in picobarns) but without 
any additional factors such as flux-factor l/(2s'), etc. In the following subsections we will 
explain how to use the interface and describe in more detail the external libraries supplied 
with KoralW. 

3.2.1 Interface 

A predefined interface, now included in KoralW, will activate user-supplied routines with 
the help of the Key4f key. For Key4f =0 no external matrix element is included and for 
Key4f =1 it is active. The new position of the weight switch KeyWgt is also introduced. 
For KeyWgt=2, the generation is done in two steps: in the first step the program generates 
const ant- weight events according to CC03 distribution and in the second step the external 
weights are calculated and transmitted to the common block /wgtall/ without rejection, 
i.e. the user gets variable-weight events. 

In order to avoid possible overwriting errors in the execution, the communication of 
the libraries with the main generator is done through dedicated interface routines. On 
the side of the generator the file amp4f_ini.f in the directory interfaces contains the 
"buffer" subprograms. On the user side, his/her own directory has to replace the directory 
ampli4f . The following two routines have to be provided by the user: 

1. ampini(xpar,npar), which initializes the external matrix element library. The 
standard input parameter array xpar can be used there for initialization purposes 
(entries above 4000). The additional parameters are made available to the user by 
calling the subroutine masow(sin2w,gpicb,ainaf ) with sin2w being the sin^6'i^, 
gpicob the conversion factor to picobarns and amaf (20) the matrix of fermion 
masses. 

2. amp4f (ql , if bml , q2 , if bm2 , pi , if 11 , p2 , if 12 , p3 , if 13 , p4 , if 14 , wtnid4f , wt4f ) 
should calculate the new matrix element squared wtmd4f, fully normalized, but 
without any additional factors such as e.g. flux factor l/(2s'). The arguments 

ql , if bml , q2 , if bm2 , pi , if 1 1 , p2 , if 12 , p3 , if 13 , p4 , if 14 denote respectively four- 
momenta and identifiers (according to the PDG conventions |2^) of the initial-state 
effective beams and the final-state fermions before the final-state bremsstrahlung 
generation. In fact the four-momenta of the final-state fermions are always supplied 



^ It is important to stress that this practical factorizabihty of the algorithm is due to the approximation 
in the treatment of complete bremsstrahlung. Omitted non-factorizable corrections may produce severe 



effects in certain final-state configurations, see Section 3.3.2 for details. 
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by KoralW to amp4f in the same order as in the /cms_ef f jnomdec/ common block. 
The additional vector of weights wt4f (i) , i=l , 9 may optionally be filled in by the 
routine amp4f . It is not used in the program but only transmitted to the KoralW 
optional weights common block /wgtall/ as wtset(40+i). The wtmd4f is put into 
wtset(40). 

3.2.2 All Four-Fermion Library of GRACE v. 2 

A working example of the above interfacing philosophy is the library of all four-fermion 
matrix elements. It is produced by the package for automated calculations GRACE 



version 2 [|12[ and provided for KoralW by courtesy of the Minami-Tateya Group of KEK. 
It consists of a set of fully massive matrix elements for all possible four-fermion processes 
in the e~^e~ collisions. Let us describe some features of the package: 

1. The CKM matrix is diagonal. For the non-diagonal CC final states the internal 
CC03 massless matrix element is used as a temporary fix. For the non-diagonal, 
doubly-CKM-suppressed MIX-type final states, i.e. simultaneously CC- and NC- 
type {ussu, ubbu, cddc, cbbc), the external matrix element (it is without CC03 graphs 
due to the diagonal structure of the CKM matrix) is incoherently added to the 
internal one (containing CC03 graphs only) and the interference of these two is 
neglected. 

2. It can be downgraded to the CC03 case with the help of the dipswitch iswitch=0 
in the amp4f routine. 

3. There are Higgs exchange-graphs included, activated by the dipswitch jhiggs in 
the file grc4f _init/setmas_koralw.f , cf. GRACE manual |T^ for details of the 
model. 

4. The gluon-exchange graphs can be activated with the help of the dipswitch jgluon 
in the file grc4f _init/setmas_koralw.f , cf. GRACE manual |T^ for details. 

3.2.3 CC-All Library of GRACE v. 1 

There is yet another external library distributed with KoralW. It is a library of massive 
matrix elements for the CC-type processes. It is based on the earlier version 1 of the 
GRACE package [0. The important difference, as compared with the previous library. 



is that complete set of scripts can be found in here to generate the library from scratch, 
i.e. to build it by the algebraic program GRACE and then custom-fit for the KoralW. 
It requires the GRACE v. 1 package to be installed. Also, the complicated structure 
of makefiles and scripts may pose some portability problems. We have developed the 
package on an HP735 computer running HP-UX v. 9.x. 

Currently, this library is obviously inferior to the GRACE v. 2 library, and we do 
not recommend using it for the standard applications. On the other hand, for an ad- 
vanced user, it gives a possibility to define his/her own modified "Standard Model" at the 
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Lagrangian level (within the GRACE system) and then directly have a ready-to-plug-in 
code generated. For this and some other reasons we decided to keep it in the distribution 
version. We recommend the interested users to contact us directly for more information 
on its use. 

3.2.4 Cuts, Instabilities and Precision 

The following set of primary cuts at the low level of MC generation is introduced in 
KoralW v. 1.42 (in subroutine selecto in korww/karludw.f): 

1. on the invariant-mass squared of produced pairs in final states with two or more 
electrons/positrons, e~^{pi)e~{p2)f{p3)f{p4,): (^1+^2)^ > arbitrl, (^3+^4)^ > 
arbitrl; in the case of e^ {pi)e~ {p2)e~^ {p3)e~ (pi) final state, additionally, {pi + 
P4,Y > arbitrl, {p2 +^3)^ > arbitrl are required; 

2. on angles of charged fermions with respect to the beams: 
Z{pi,hearaj) > themin [rad] 

3. on transverse momenta of visible charged fermions p'^ > arbitr with p'^ defined as 
Py = Yli (pf(l) +P?(2)) summed over all charged fermions with the angles greater 
than themin with respect to the beams; 

4. on transverse momenta of photons for e+e~//-type final states: 
(E-P.(l))'+(E-P.(2))'<arbitr2 

The default and recommended values of these cuts are: 

arbitr = 600 GeV^ for minimal visible p"^, 

arbitrl = 8 GeV^ for minimal invariant mass in e'^e" ff, 

themin = lxlO~^ rad for minimal angle with beam, 

arbitr2 = 300 GeV^ for maximal p^ of photons in e^e~//. 

Let us now turn to the basic question: Do we need these cuts and why? Before we try to 

answer it, let us stress that these cuts can be removed completely. Both the presamplers 

cover the entire four-fermion phase space, no part of it is cut-out during the generation 

process. The above cuts are imposed after the generation is completed, but before the 

matrix element is calculated. Now, coming back to the basic question, there are two 

reasons for non-zero recommended values of the cuts: 

• Technical problems in the evaluation of the external matrix element. It can be 
numerically unstable for some extreme four-momenta configurations (when some of 
the invariants become very small). 

1. For the CC-type processes we have not encountered any instabilities in any 
region of the phase space while using presampler No. 1 in the entire phase- 
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space, so no cuts are needed^. In the case of presampler No. 2 we noticed 
some residual instabilities that can be cut out with the angular cut themin 
= 1x10^^ rad. This cut influences the total cross-section at the level of 1-2 
per mille and therefore can be neglected within current precision requirements. 
(The presence of instabilities only in presampler No. 2 is due to the fact that 
this presampler is more oriented on the singular configurations and in a sense 
reaches "closer" to the singularities than the other one.) 

2. The other two cuts are for the processes classified as NC-type. The instabilities 
show up either for small transfers in the t-channel-dominated configurations 
(cut out by arbitr) or for small invariant masses of fermion pairs (cut out by 
arbitrl). Note that the arbitr cut, although defined by transverse momenta, 
is intended to assure that the energy of the (other and undetected) final-state 
electron, even if lost in the beam pipe, is noticeably smaller than the beam 
energy. 

3. The above three cuts have proved to us to be sufficient to deal with the Born- 
level instabilities in the matrix element (see previous footnote). Unfortunately 
the picture gets more complicated in the presence of ISR bremsstrahlung. 
Namely, with the emission of the transverse photons the effective CMS frame 
that we use to define the four-momenta for the matrix element calculation gets 
rotated with respect to the LAB frame in which the cuts are imposed (see 
Ref. |l[] for details on the construction of this effective frame). As a result, 
a "large- angle-event" in the LAB frame can lead to the highly coUinear and 
numerically unstable matrix-element calculation in the effective CMS frame. 
For that reason the arbitr2 cut is introduced along with the relatively high 
values of arbitr and arbitrl cuts. 

In principle, the instabilities of points 1 and 2 can be cured with the help of 
quadruple-precision arithmetics. It is an easy task on the IBM or SGI worksta- 
tions as they support C0MPLEX*32 type and the upgrade can be done easily with the 
help of the compiler options. This cures the technical side of the problem. How- 
ever, the problem described in point 3 will remain unsolved, as the cross-section 
(and weights!), although numerically stable, will be huge and very likely physically 
wrong in the presence of transverse photons. It is due to the ill-defined simulation 
of the missing complete 0{a) matrix element and the use of the effective beams 
technique. 

Physical problems. At the end of the last point we already presented one aspect of 
physical inaccuracy due to the lack of a complete 0{a) matrix element. The other 
related fundamental problem is due to the fact that the bremsstrahlung process is 



^ Of course we cannot guarantee that these instabihties would not show up with sonic combination of 
input parameters, with sufhcicntly long series of generated events or even on some platforms other than 
the ones we tried. 
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treated in the Initial State Radiation approximation. It means that only the s- 
channel type radiation is generated. In the case of electrons in the final state, one 
often encounters for example the Bhabha-like configurations that are dominated by 
t-channel-type diagrams. The radiation, however, is neither generated according to 
nor corrected for these configurations. This can lead to severe inaccuracies in the 



results. We refer the reader to Section B.3.21 for more information. 



To summarize, the technical instabilities in the Born cross section can be cured either 
by going to the quadruple precision or by implementing relatively modest cuts, weaker 
than the recommended ones. The inclusion of bremsstrahlung makes the picture compli- 
cated. Because of the limitations of the effective beams technique, even the recommended 
strong cuts do not remove the instabilities completely. On the other hand, even if the 
quadruple precision is used one is still limited, in loosening cuts, by the physical meaning- 
fulness of the results, and, on the technical side, by large overweights due to overestimated 



cross-sections. This subject requires further study, see also Section |3.3.2 . 

Let us also note that in the case of the CC03-type presampler all the cuts are reset by 
the program to no-cuts values. 

Together with the predefined cuts in the selecto routine, we provide another routine, 
user_selecto, in which the user can implement his/her own cuts. This routine is located 
in the same file korww/karludw . f . If the user requires some stringent additional (pre) cuts 
they can be placed there. This may increase the speed of the program (depending on how 
many events are rejected by the cuts) as the matrix element will not be calculated for 
unaccepted events. The user interested in the complete or nearly complete phase-space 
coverage should start with setting all cuts to zero and inspecting overweighted events. 
(The arbitr2, if set to a negative or zero value, will be reset by the program to its 
maximal value, i.e. the value of the s parameter). The maximal weights for rejection can 
be adjusted in the input parameters later on. 

As an example of such a routine, we provide the routine user_selecto_canonical 
with the canonical cuts as defined at the LEP2 Workshop ^. With its help we have for 
example reproduced, within an accuracy of a few per mille, the results given in Tables 6-8 
(see p. 241 of Vol. 1) in Ref. Q (with the exception of entry 14 of Table 7). 

3.3 Bremsstrahlung 

In the present version of our program we upgraded the ISR exponentiated matrix element 
with third-order LL corrections. We know that these effects are small, especially for the 



Yennie-Frautschi-Suura-type exponentiation we use, see Ref. |26]; this new contribution 
may nevertheless be useful for estimating theoretical systematics due to QED higher 
orders. The reader should keep in mind that the new third-order LL corrections influence 
longitudinal momenta of ISR photons and add nothing new in the pt distributions of the 
ISR photons. Let us write the master equation for the total cross section, with new terms 
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due to third-order LL corrections: 



a 



|:;^/nf^(nf^5ta,p.,*.)')««(p,+p.-t,.-|:*.)e 



exp Ua^B + 






sip,,p2,k)ii-ei 






_^yP?ip^,qf,h,k,) 



i>j 



b [ki)b[kj ^ 



sr-^ P3 [Pr i Qs ) ^»? kj, ki^ 

^^^ ~s{K)~s{k,ys{h) 



where S{pi,p2,k) = —{a/4:n'^)(^{pi/kpi) — {p2/kp2)) is the real photon infrared (IR) 

n 

factor and Qf^ = IT ^ (^f " ^V^/^) cuts out the singular IR region, already included to 

all orders in the YFS form factor ||l],0. In Eq. (Q) we use essentially the same notation 
as in Ref. ^, except for p^ = TZpr, r = 1,2 and qf- = TZqs, s = 1, ..., 4. Note that the 
"reduction procedure" 7^ does not act on photons. The new third-order non-IR functions 
Pi , / = 0, ..., 3 are defined below. In our 0{a^) LL calculation all /5's are proportional 
to the basic lowest-order (Born) distribution 
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in which fermion momenta are adjusted (with 7?.-procedure) in such a way that the Born 
matrix element A^Bom is calculated at the reduced centre-of-mass energy s' = (pi + P2 — 

3.3.1 ISR up to Third Order 

The 0(0;'") functions (3^ , see 0, are residuals of the removal of IR virtual and real 
singularities. They are therefore IR-finite. In the YFS scheme, they are obtained from 
the (9(a'') "raw" differential distributions, which originate from the Feynman diagrams. 
This is also the case in the 0{a^) YFS-exponentiated matrix element of BHLUMI |^^ or 
YFS2 1^ and YFS3 |^. Sometimes the more "economical" source of raw differential 
distributions is not the Feynman rules but the LL approximation, especially for 0(7^) 
and 0{'j^) corrections, where 



7 



a 



vr 



In- 



1 



mi 



(6) 



Here, in our calculation, we employ the LL approximation up to 0(7^). 

The first step in our procedure of deducing the third-order LL exclusive distribution is 
to examine triple convolution of the Altarelli-Parisi (AP) non-singlet evolution equation. 
This is done separately for two e^ incoming lines and the convolution of the two results 



1^ 



is examined once again (truncating the final result up to 0{'-f^)). It can also be done 
in one step, using a triple convolution of the AP equation, with a "double" AP kernel 
acting on both e^ lines. All the above is a one-dimensional exercise, i.e. variables in the 
AP equation are interpreted as ratios of the photon energies to the e^ energies, and the 
photons are assumed to have exactly zero transverse momentum. The above calculation is 
rather straightforward, quite similar to that done in Ref . pO| , and we shall here skip this 



part. With the explicit longitudinal momentum of up to 3 photons, we supplement them 
with the smooth pt distribution, getting in this way what we call the "LL ansatz" for 
full differential distributions, with up to three real photons. In the ansatz we are careful 
to reproduce all correct soft-photon limits, i.e. for any configuration with any number of 
photons being soft. In the above construction, virtual up to 0(7^) corrections are kept 
all the way through the construction procedure, and the cancellation of IR singularities 
is kept perfect. The resulting (9(7^)-finite pt ansatz for n^ = 0, 1, 2, 3 photons reads as 
follows: 

^fo] = {1 + 25(1) + l[2B{l)r + l[2B{l)f^ b^, 

D^^ih) = 5(A;i)x(«i, /?i){l + ^[35(1) + B{z,)]' 

+ i [(25(1))^ + (2i?(l))(i?(l) + B{z,)) + {2B{z,)r] 



''0 ' 



D^^{h,k2) = Sih)Sik2)x{ai,(3,)x{a;,(3;) |^ + ^[AB{1) + B{z,) + B{z,Z2)]\ b^, 

D§^{h,k2,h) = Sik,)Sik2)Sik,) xiauPi)xK,(3*2)xK,P;) &?, 
;^(a,6)^((l-af + (l-6f)/2, 

B{z) = ^A{z), A{z) = ^ + \ne-\nz, 

(7) 
where Oj and /3j are the familiar Sudakov variables 

tti = ihPi) / iPiP2) , Pi = {kiP2) / {P1P2) (8) 

while "starred" variables take into account the effects of energy loss due to prior emissions 

«; = «2/(l-ai), f3*2=f32/{l-Pi), .g. 

ag = 03/(1 - ai - 02), Ps = /33/(l - Pi- P2), 

and we define z-variables as 

^1 = (1 - ai)(l - A), z,Z2 = (!-«!- a2)(l - A - P2). (10) 

The e variable is an infrared cut-off (an IR-regulator in the AP kernel), which defines the 
lower limit of photon energy in our phase space, k^ > eA/i/2. Note also that the fact that 
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7 comes from the angular integration over S{k) is essential for the exact IR cancellations 
between the virtual and real photons. The first step on the way to (3}J is to remove IR 
virtual corrections (because they are already present in the YFS form factor) 
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and the resulting distributions (programmed this way in our code) read as follows: 
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We assume tacitly that Bose symmetrization is done in the above expressions, for n = 2,3. 
Finally, the real IR-parts are subtracted properly as follows: 
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- S{ki)S{k2)S{k3)P^''\ 

The above subtraction is done numerically in the program. This completes the definition 
of our 0(7^) matrix element. Let us note finally that the above parametrization was 
essentially copied from the %% Monte Carlo for fermion-pair production ^1 . 



3.3.2 Non-ISR Corrections 



As was already mentioned in Subsection |3.2.4| , in the processes with at least one e"^e~ 
pair in the final state some problems are encountered when we try to include the leading 
QED corrections by convoluting the s-channel-type ISR bremsstrahlung with the Born- 
like e^e~ —^ e^e^ ff process, calculated in the "effective beams" rest frame, i.e. in the 
rest frame of the incoming e~^e~, with reduced energies after emission of the ISR photons. 
Generally, such an approach leads to two kinds of problem^: (1) too much radiation 

^ Here, we do not discuss the problem that arises when any of the electrons/positrons are lost in the 
beam pipe; this requires special treatment. 
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is generated, particularly in the high photon-py range, and (2) the event weights are 
submitted to huge fluctuations. All this is because the cross section in the above 4-ferniion 
production channels is dominated by the low-angle Bhabha-like processes, where the t- 
channel 7 exchange plays a crucial role. As a result, the photon radiation is governed by 
ln(|t|/mg) rather than by ln{s/ml). These two "big logs" differ by a factor ln(s/|t|), which 
can be sizeable at low angles (in fact, it can itself be regarded as a "big log"); this factor 
can explain an excess of the ISR bremsstrahlung when generated as in the s-channel. The 
reason for too many high-p-r ISR photons is that in the low-angle Bhabha-like process 
the interferences between the electron and positron lines become very small, while the 
destructive interference between initial- and final-state emissions becomes strong, see e.g. 



Ref p^. As a result, the emission of high-py photons is strongly suppressed, while 
in the current implementation of the s-channel-type of the ISR photon generation this 
interference is absent and cannot damp the pr of photons. The ultimate solution to this 
problem is the inclusion of the bremsstrahlung generation of both the s-channel and t- 
channel in the MC algorithm. This, however, is not easy to implement into the current 
structure of the program^. 

The second kind of problems, i.e. the strong fluctuations of the event weights, are 
partly connected to the first one. Namely, the emission of the high-p-r ISR photons can 
boost the fermions to very small angles configuration, where the "effective" Born-like cross 
section is huge, even if, in the LAB frame, the final electrons/positrons are emitted at 
moderate angles. This can also cause some numerical instabilities in the matrix element 



calculation, as has been discussed in Subsection |3.2.4 This deficiency can be partly 



cured by solving the first problem, i.e. by damping the high-pT photon radiation. But 
even after this is done, some huge weight can still appear as some residual high-pT photons 
still remain, e.g. in the configurations where only one of the electrons/positrons is emitted 
at low angles. 

We could try to overcome the above difficulties by imposing some additional cuts on 
the e~^e~ ff events, as described in Subsection |3.2.4| . Our experience shows, however, that 



even the presence of strong cuts does not solve the problems completely, although it can 
stabilize the event generation considerably. Besides, some of those cuts may not be easy 
to implement in the real experiment. 

In the mean time we are preparing an intermediate solution, without going into the 
t-channel bremsstrahlung generation and/or the exact (off-shell) 0{a) matrix element, 
which eliminates extra cuts. We are constructing an additional weight that would emulate 
the main features of the t-channel bremsstrahlung. As a guideline, we use the YFS 
approach to the radiative Bhabha-like process, looking at the behaviour of the hard 
bremsstrahlung matrix element at low angles. Since this solution is not yet ready, we 
do not include it in the current version of KoralW. Instead, we add a dummy routine 
eexx_wt_cor (the file model/eexxcor .f) as an "open slot" where the real one will be 
"plugged in" when completed. 



®We consider this solution for the future versions of KoralW. 
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3.4 Additional Corrections (CC03-Based) 

In this section we will describe the additional corrections that are imposed on top of 
the Born-level WW matrix element. All these corrections have one common feature - 
they are justified/calculated for the CC03 subset of graphs only. However, it is of great 
importance to have them incorporated, all at the same time, into the complete charged 
current (CC-all) matrix element in one way or another. A detailed study on this subject 



has been presented in Ref. [11] on the example of the KoralW and grc4f codes. On the 
one hand, it proved that the "common-sense interpolation" does make sense, provided 
one performs some minimum set of cross-checks. On the other hand, it showed that 
indeed it is a "code-dependent" procedure and the differences can be of the order of a few 
per mille. The possible choices for adding corrections are for example: additive versus 
multiplicative method or the amplitude level versus the squared amplitude level. In the 
following we will describe some of the effects: naive QCD, the non-diagonal CKM matrix, 
the Coulomb correction and the anomalous couplings. In the future we hope that this 
list will be extended by the complete 0{a) corrections for the on-shell WW production, 
available in the literature [p2-R6[ and already implemented by us in another Monte Carlo 



code - YFSWW |37, 38 



3.4.1 Naive QCD and CKM Matrix 

For the lyiy-type final states the naive QCD correction (NQCD) and the non-diagonal 
CKM matrix are applied as a multiplicative correction to the whole matrix element. In 
principle these corrections are justified for the CC03 subset of graphs only. However, 
as the background graphs contribute very little at LEP energies, we follow the standard 
approach and apply them to the entire CC-all matrix element 

\Mccau{^,J)\' = \M^cSZr%j)\' ^%QCDi^) ^%QCDij) (14) 

where \Mccaii{hJ)\'^ is the matrix element for the {i,j) decay channel of the VT-pair as 
used by KoralW for calculating weights and |iVf^'^^™"'(i, j)p the external matrix element 
supplied to KoralW by the external library for the {i,j) decay channel of the VT-pair. 
We note in passing that in the CC03 case, the CKM and NQCD corrections are already 
included in the iMccosih j)\'^ and there is no need for the correcting factor ^%qcd i^ 
this case. The correction itself is calculated for each W as the ratio of appropriate W 
branching ratios: 

rvCKM (■. _ Br(i) Brojeu,) 

with Btq being the "bare" W branching ratios (1/3 and 1/9). The physical W branching 
ratios are set either directly to fixed values (KeyBra=l), corresponding to 0:5 = 0.12, or 
calculated according to the formula of P^ based on the CKM matrix (KeyBra=2): 



Br{q) = \\VcKM{q)?^^^. Br{l) = \^^. (16) 

3 TT 3 TT 
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In either case the Tw is forced to be recalculated by the program. The third option 
KeyBra=0 is equivalent to KeyBra=2 with the diagonal CKM matrix (1/3 for quarks and 
1/9 for leptons) and F^y not recalculated. The formula used for T]y recalculation (default 
if the input value of T^r is negative) is the following: 



^ ^^^^^^^ (17) 

Note that, unlike the situation in the previous versions of KoralW, here, the W branch- 
ing ratios are used solely for the normalization of the matrix element. 

For the ZZ-tjpe final states, NQCD is applied directly as a multiplicative correction 

\M^Call{^,J)\' = \M^^%lT\^,J)\'^NQCD{^)^NQCD{J), (18) 



^NQCD{q) = l + —, ^nqcd{1) = 1, (19) 

vr 

to the whole matrix element, for each generated quark pair. The Z-width is however 
not corrected automatically, and the user must take care of it by himself/herself (since it 
enters the program through the input parameters, it can easily be set to the appropriate 
value) . 

3.4.2 Coulomb Correction 

The Coulomb correction is taken from Ref. [Q, Eq. (9). It is the first-order formula: 



1 Q — 4 /If^ 

/ = Us' - 2s{s, + s,) + {s, - s,n E = "—iP^, (21) 

4s 4:Mw 




E^ + TI^-e] - iJ]^Mw Ue^+TI, + e\ (22) 



In the case of the CC03 matrix element, "^Icoui is a simple over-all multiplicative correction. 
In the case of the CC-all matrix element, it is implemented in the code in the form of an 
additive correction based on the CC03 matrix element only, and Eq. (p^ ) now becomes: 



\McCau{l,j)? = \Mh%Zir\^,j)?'^%CD{m%QCD{j) + WcCoS, j)\\'^Coul ' I) ■ (23) 

An identical correction is applied to the CC03 semi-analytical formula. 
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3.4.3 Anomalous Couplings 

In the case of the complete charged current, CC-all, matrix element, the anomalous 
couplings are implemented in the code in the form of an additive correction. Equations 
(|l^ and (^) take on their final form: 

+ {\M^So3ihJ)ncoul-\Mcco3{^,J)\'), (24) 

which gives the complete description of the matrix element of KoralW, including all the 
"additional" effects. 

In the current version of the program we implement three parametrizations of the 
anomalous WWV couplings {V = 'j or Z). A particular parametrization can be chosen 
by the user with the help of the input parameter switch KeyAcc. Values of these couplings 
are transferred to the program through the input parameters vector xpar, and they have 
to be set up appropriately by the user in the initialization mode. 

The parametrizations of the triple gauge boson couplings are the following: 

• KeyAcc=0: the Standard Model (SM) (non-anomalous) couplings. 

• KeyAcc=l: {gf , Ky, Xv, 9^ , dhf^v^K} ^ 'C, (^ = 7,^) 

The above set represents the most general parametrization of the WWV couplings 
that can be observable in the process where the vector bosons couple to effectively 
massless fermions |^D|; see also Ref. |^T| and references therein. In general, there 



are 7 complex-number-valued couplings for each 7 and Z, so altogether one needs 
to supply 14 complex (or 28 real) numbers in this case. This can be done through 
the xpar vector entries: xpar (21-57), see Table for more details. The SM values 
for these couplings are: 



V 

9i = i^v 



V V ~ \ n 

9i =95 = nv = Xv = 0. 



This parametrization has been present in KoralW since the version 1.03 and is back- 
ward compatible. The new parametrizations added in this version of the program, 
for KeyAcc=2,3, are described below. 

KeyAcc=2: {Sz, x^, Xz, y^, yz} e M 

This set of 5 real-number parameters represents deviations of the C- and P-conserving 



couplings from their SM values; see Ref. [41| and references therein. It can be related 



to the previous parametrization as follows: 



9j = 


= 1, 


9l- 


= 1 + tan6'vy (5z, 


r\jiy — 


= 1 + x^. 


Kz -- 


= 1 + tan6'vy(a;z + bz 
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^7 — y-y, 




Az = yz, 




\v = 9l- 


-9\ 



Ky = Ay = 0. 

The values of the above 5 parameters are sent to the program through the xpar 
vector entries: xpar (61-65); see Table ^ for details. For the SM they are all zero. 
This parametrization is convenient for studies of the C- and P-conserving WWV 
couplings. 

KeyAcc=3: {avy?i, «b<^, «H/} ^ M 

This parametrization of 3 real-number- valued couplings corresponds to a direct ex- 
tension of the SM formalism in terms of a linear realization of the symmetry, which 



can be achieved if a relatively light Higgs boson is assumed to exist, see Ref. |^T| and 
references therein. Its relation to the most general parametrization (for KeyAcc=l) 
reads: 



9l- 


= 1, 


9l 




rVj'y 


= 1 + awif, + dBcf. 


i^z 


s2 

= 1 + aw<f, 2 

Cy/ 


A7 


= \z = OL-Wi 


Ay 


V V ~ 

= 9i =95 = i^v 



Ay = 0, 

where sw = sin6'vy, cw = cos 6w- 

Values of the above 3 parameters are sent to the program through the xpar vector 

entries: xpar (71-73), see Table |]for details. For the SM they all are zero. 

3.5 Renormalization Schemes 

The choice of renormalization scheme is a direct result of the freedom in the renormal- 
ization procedure of the electroweak sector. One can choose a different set of input pa- 
rameters used to express the cross-section. There is a large number of different schemes 
available in the literature; see, for instance [Q. All these schemes are more or less edu- 
cated, but ad hoc inclusions of some EW radiative corrections, and are done by hand on 
the "bare Born" by means of improving its parameters. Therefore they can in principle 
vary from the "bare" (not corrected) result and, among themselves, as much as 15%; see 
for example [Q. 

In KoralW we have implemented three different schemes. Any of them can be chosen 
with the help of the input parameter KeyMix: 

• KeyMix=0 "LEP2 Workshop scheme" : 

a^y = (input), sin^ 6'ty = 7ravy/(y2G^M^), gf^ = 47rQ;vy/sin^6'vy; 
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• KeyMix=l "G^ scheme": 

aw = V2Gf,M^ sin^ Ow/t^, sin^ Ow = I - M^/M'^, g"^ = Aixaw/ sin^ Ow] 

• KeyMix=2 "bare Born" (for tests only): 

o^w = otQED = 1/137, sm^ 9w = 1 — M^/M|, g"^ = Attuw / ^"vn^ Ow 

Note that this key has changed its meaning with respect to the earher versions of 
the KoralW code, where it simply changed the definition of sin^ 6^, leaving changes of 
the other input parameters directly to the user. Now it does the complete redefinition 
of the required input parameters according to the chosen scheme. The reader may have 
noticed that the options KeyMix=l and KeyMix=2 have changed since the previous version, 
whereas the setting KeyMix=0 remained unchanged. 

Finally, let us comment on the issue of the recommended setting. The ultimate cri- 
terion of choosing the scheme would be the comparison with the exact electroweak cor- 
rections. This has been studied only for the on-shell case 0, and no unique answer was 
given. Favoured were, however, the "Workshop" or "G^" type schemes. In fact, it would 
most likely be impossible to distinguish them experimentally, as for example these two 
favoured schemes differ at the per mille level, far less than the correction itself. 

3.6 Colour (Re) Connection 

3.6.1 Colour Connection 

The final state of KoralW generation consists of four fermions and an arbitrary number of 
real photons. Later, external libraries are called, if necessary, to perform, for instance, the 
r decay and/or hadronization of the quark pairs. In most cases, the information stored 



by KoralW in the standard event record COMMON /HEPEVT/ |g5[ is sufficient. 

However, there exists a sub-class of four-quark final states where the intermediate state 
can consist of two distinct configurations of colourless quark-pairs, e.g. for the uu, dd final 
state the colourless pairs of quarks chosen for hadronization can be either ud and du or uu 
and dd. This ambiguity has to be determined by KoralW, so that appropriate information 



can be passed to the hadronization package JETSET [15 



In our program, the random choice is made with the help of the weight 

|Mi|2 



wcc 



2 I I A/f.12 



|Mi|2+|M 



where Mi and M2 represent spin amplitudes squared (and appropriately summed over 
spin degrees of freedom) for the two possible colour-singlet configurations. The choice is 
performed with the help of a call to the routine spdetx(ireco). The output parameter 
ireco=0 , 1 denote that the colourless object should be formed either from the first-second 
and third-fourth quark pairs in the final state (as located in the COMMON /HEPEVT/) or 
from the first-fourth and second-third quark pairs, respectively. 

As an input for the routine spdetx(ireco), the appropriate spin amplitudes hidden 
in internal common blocks of the GRACE routines are used. The routine spdetx(ireco) 
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consists of the routine spdetc generated by the GRACE package [^, which is adapted 
to our purposes. 

3.6.2 Colour Reconnection 

Since in the LEP2 energy range the average distance between the W~^ and W~ decay 
vertices is much smaller than the typical hadronic size, the fragmentation of two W^s may 
not be independent. One of the physical effects related to this problem is the colour re- 
connection between decay products of different W^sf\. Several models have been proposed 
to describe this phenomenon and estimate its influence on the 1^-mass measurement; see 
e.g. Ref. P] and references therein. The current version of KoralW does not include any 
of these models; it uses instead the user-supplied value of the colour-reconnection prob- 
ability PReco (see Table P to generate quark-antiquark colour objects from the decay 
products of two different W^s. These objects are then processed by the JETSET routines 
to result in the final-state hadrons. 

This, rather simplistic, approach to the colour reconnection problem does not deal, 
of course, with all the aspects of this phenomenon, but it can be used for some simple 
estimate of its influence on physical observables. In detailed studies, however, this has to 
be confronted with the results of dedicated models. 

3.7 Bose-Einstein Effect in Hadronization 

It is generally expected that in the double hadronic decay of the H^-pair, because of the 
space-time overlap of the two hadronization processes, the two W's cannot be treated as 
completely independent. Consequently, one may see experimentally some effects due to a 
"cross-talk" between the decays of the two W's. One of the possible effects (in addition 
to colour reconnection) could be a deformation of the hadron distributions due to the "co- 
herence/interference" effects in the hadronization process called the Bose-Einstein effect. 
This effect is of interest in itself, on the other hand it also can obscure the measurement 
of the W mass in the decay of PF-pair using four-jet final states. There are a variety of 
models, see 0, describing the BE effect. In particular, the JETSET package implements 
one of these models with the help of routine LuBoEi. The LuBoEi routine introduces 
additional smearing of the momenta of the hadrons, after standard hadronization, so that 
the Bose-Einstein effect is reproduced. In KoralW we include an example of the alter- 
native model, see Ref. [^ , where this additional smearing is done not by means of direct 
manipulation of the hadron momenta, but with the help of an additional special weight. 
In principle, this method is safer because it does not introduce spurious long-range cor- 
relations of hadrons; it is therefore more relevant to ly-mass measurement from jet-jet 



^ In the general case, not restricted just to the W-pair production and decay at the CC03 approxi- 
mation, the quark pairs, as explained above, can originate also from decays of Z's or virtual 7's, and the 
quark-pair colour-singlet ambiguity can already exist in that step. Independently, in the second step, the 
colour arrangement may need to be redefined, i.e. reconnected, owing to the final-state giuonic/hadronic 
interactions. 
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effective masses than the LuBoEi procedure. The BE weight itself is programmed in the 
C++ class BEwtMaker in subdirectory B.E./src. (It uses the additional C++ functions 
hidden in the file partitl.C.) The code is rather compact, the bulk of the code is the 
construction of "clusters" of pions of the same sign, which are in some sense close to one 
another in the phase space (their relative hyper- velocity is below certain maximum value). 
The subdirectories B.E./src and B.E./fig contain the complete tool-box for analysing 
the BE weight and finding out how big the BE effect is in the fitted W mass. For more 
details on this analysis of the BE effect programmed in the subdirectory B . E . , see Ref. 0. 

3.8 Energy Distributions 

In this section we describe part of the code responsible for the generation of the s' photonic 
variable and decay channel type. The important observation is that the ISR photons and 
final-state fermions are coupled solely by the s' variable. This variable, on the one hand, 
determines the amount of energy radiated by the photons and at the same time provides 
the effective CMS energy of the final fermions. Its generation is therefore an important 
step in the algorithm. 

As usual in the Monte Carlo algorithms, one can start from anything like the s' 
distribution, at the end of the MC generation procedure, this dummy distribution will be 
transformed into what we want and the original one will get eliminated (by reweighting 
and/or rejecting MC events). The price for starting with a bad distribution is, as usual, 
in the effectiveness of the program. Since the cross-section varies by orders of magnitude 
with the energy and between various final states, the overhead in efficiency can be really 
big. For that reason we decided to pretabulate s' distributions. It is done for each final 



state separately and for the "standard" set of cuts, as given in Sect. |3.2.4 . 

We see the potential disadvantage of this organization as well. In the case of cuts dif- 
ferent from the "standard" ones, pretabulated spectra would gradually become inefficient. 
In our opinion, the gains due to the precise tune- up outweigh, however, the loses. In fact, 
one can always generate another set of files with new spectra, if one is interested in very 
different cuts. We have created semi-automated tools for this task (not distributed with 
KoralW 1.42) and encourage the interested user to contact us directly. 

There is yet another trick that we used here. It seems natural that the pretabulated 
distributions are Born cross sections (as a function of CMS energy). For the s' generation, 
they will be convoluted with the structure-function-type photonic density, with proper soft 
singularity, and they are used directly for flavour generation. However, the complicated 
four-fermion phase space leads, in some places, to events that are over-weighted with 
respect to the global rejection weight. A global increase of the rejection weight would 
cause an unnecessary loss of effectiveness in the case of a run with constant weight. 
The pretabulation procedure distributions provide, as a by-product, an efficient tool for 
compensating the above effect: with their help, we may provide the system of additional 
compensating weights, depending on the flavours and the CMS energy. In order to do it 
in a transparent way, we therefore pretabulated the maximal weights, in addition to the 
total cross sections. And in fact it is these maximal weights that the program uses for 
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the generation of s' at the crude level in the CC/NC-all mode. In the CC03 mode, the 
total cross section is used and the weights are well behaved. 

We provide also an option for generating the distributions from the user-supphed 
function. The function is called phot_spec_crud and the appropriate dip-switch i_f ile=0 
is located in karlud routine. It is capable of reproducing the distributions used in the 
previous versions of KoralW. 

At the technical level the organization is the following. The pretabulated files, located 
in the directory data_files, are: data_wtmax.f it . smpl, data_wtmax.f it . smp2 and 
data_xsect .f it. They contain distributions of the maximal weight for two presamplers 
and the total cross section, respectively. They are tabulated as a function of the total 
energy for each final state separately. The tabulation is done in the energy range from 50 
GeV to 250 GeV; beyond these limits, a flat-distribution-matching value at 250 GeV is 
used. The give_phot_spec_crud routine reads these files and stores the data in some local 
internal variables. KoralW uses another set of matrices with different binning (currently 
the finer one): prob_*_* (here * denotes a wildcard for a part of a variable name), and 
a different energy range, from some emin set to 1 GeV in the korww routine, to the 
emax=A/s. The arrays prob_*_* contain channel- dependent distributions for the channel 
choice as well as an inclusive distribution for the s' generation. 

3.9 Semi-analytical Distributions 

The package for semi-analytical calculations distributed with KoralW integrates the CC03 



matrix element, based on the formulas of Ref. |^. The calculation is done in massless 
approximation. Additional corrections due to NQCD, the non-diagonal CKM matrix and 
the Coulomb effect are included in the formulas. It is done exactly as in the Monte Carlo 
generator, so the two nicely cross-check each other. The Initial State Radiation is also 
included. It is done in the LL approximation, by means of the Structure Function formal- 
ism. The SFs are implemented up to the third order with the YFS-type exponentiation 



and the non- leading YFS form factor p6|,|30|. For the actual calculations we provide a 



number of different approximations of these SFs, as listed in Table |I4|. For more details 
on the integration method see Ref. [|1^ . 

In the current version of KoralW the semi-analytical part of the program is enlarged 
with two functions slwan(sl) and sls2wan(sl,s2) for the one- and two-dimensional 
distributions of the single and double W invariant masses. These functions require stan- 
dard initialization of the korwan routine. Optionally, if the KorWan input parameter 
KeyMod is increased by 10000 the calculations in KorWan are not executed and only the 
initialization is performed. 

In addition we provide two simple applications of KorWan: the average mass and 
average mass-loss calculations. The average mass is defined as (l/cr) J dvdslds2{^/sl + 
^/s2 — 2Mw)da / {dvdsids2) , with v being the photonic variable v = l—s'/s and Si, S2 being 
the invariant masses of the 'W-states" . This calculation can be performed by invoking 
KorWan with a negative value of the s-variable input parameter. Alternatively, one can 
simply use the dedicated routine mavrg with arguments as given in Table 0. The average 
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mass loss is defined as {y/s/2){l/a) J dv v da/dv with v = l—s'/s. This calculation can be 
performed by invoking KorWan with the negative KeyPho input parameter. Alternatively, 
one can simply use the dedicated routine mloss with arguments as given in Table |l^. 

4 Practical Use of the Program 

In this section we will familiarize the reader with the input and output parameters, and 
the usage of the present version of the KoralW package. We will also present two simple 
demonstration main programs using the KoralW package. Their double role is to serve 
as a useful template for the user to create his/her own main program and to help the 
user to check quickly that the KoralW generators runs correctly. We shall describe in 
detail all input parameters of KoralW. We will also give a similar information on the 
semi-analytical routine korwan. 

4.1 Principal entries of KoralW 

The principal entries of the KoralW package, which the user has to call in his/her ap- 
plication in order to generate series of the MC events, were already listed and described 
briefly in Subsection p.ll| . Here we shall add more information on their functionality. The 
calling sequence constituting a typical Monte Carlo run will look as follows: 

CALL KW_ReaDataX(' ./dataJDEFAULTS' ,1,10000, xpar) ! reading general defaults 

CALL KW_ReaDataX(' ./user . input ' ,0, 10000, xpar) ! reading user input 

CALL KW_Initialize(xpar) ! initialize generator 

DO loop=l, 10000 ! loop over MC events 

CALL KWjyiake ! generate single MC event 
EMDDO 

CALL KW_Finalize ! final book-keeping, print 

CALL KW_GetXSecMC(XSecMC,XErrMC) ! get total cross section 

In the first call of KW_ReaDataX, default data is read into the array REAL*8 xpar (10000). 
The KoralW has almost no data hidden in the source code. (This is not true for TAUOLA 
and JETSET). The file data_DEFAULTS is read first into array xpar. This file we provide 
in the distribution subdirectory data_f lies. The user should never modify it. It can be 
copied to a local directory or, better, a symbolic link should be created to the original 
data_f iles/data_DEFAULTS. The data_DEFAULTS is rather big and the user is usually 
interested only in changing some subset of these data. In the second call on KW_ReaDataX 
the user can overwrite the default data with his/her own smaller set on input data which 
are placed in the user . input file. For example the simplest input data, which defines 
only CMS energy looks as follows: 

BeginX 

*<ia>< data >< comments > 

1 190d0 CmsEne =CMS total energy [GeV] 
EndX 
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As we see, data cards start with the keyword BeginX and end with the keyword EndX. The 
comment hnes are allowed - they start with * in the first column. The data themselves are 
in a fixed format, with the address i in xpar(i) followed by the data value and trailing 
comment. The four examples of input data sets for the two demonstration programs 
KWdemo . f and KWdemo2 . f in the subdirectory demo . 14x/190gev provide useful templates 
for the typical user data. The complete set of all user data in data_DEFAULTS is described 
in very detail in Tables |l|-§|. Obviously, the user is interested in manipulating only some 
of them and will stick to default values in most of the cases. The KW_Initialize is 
invoked to initialize the generator. It reads input data from array xpar, prints them 
and sends down to various modules and auxiliary libraries. The programs have to be 
called in strictly the same order as in the above example. At this point we are ready 
to generate series of the MC events. The generation of a single event is done with the 
help of KW_Make. After the generation loop is completed, we may invoke KW_Finalize, 
which does final book-keeping, prints various pieces of information on the MC run, and 
calculates the total MC integrated cross section in picobarns. In order to obtain this cross 
section the user may call the routine KW_GetXSecMC(XSecMC,XErrMC). 

4.2 Input /Output Parameters 

As we have already explained in the previous section, the input parameters enter through 
the xpar array, being a parameter of CALL KW_Initialize. The meaning of all of them 
is given in Tables |I|-§|. 

The principal output of KoralW is the Monte Carlo event, which is just a list of 
final-state four-momenta in [GeV] units and fiavours, encoded in the standard /hepevt/ 
common block. At the present version we still provide REAL*4 version of the /hepevt/ 
common block. If the user is interested in the parton momenta before hadronization, 
then, in addition to /hepevt/ common block, they are also available through 



CALL KW_GetMomDec(pl,p2,p3,p4) 

CALL KW_GetBeains(ql,q2) 

CALL KW_GetPhotAll(NphAll,PhoAll) 



get momenta of four final fermions 

get beam momenta 

get photon multiplicity and momenta 



Alternatively, all the four-momenta are available via the internal commons /momset/ and 
/momdec/, see Tables ^ and ^ for details. 

For some special purposes, also the four-momenta in the effective CMS frame (used 
by default by the code for matrix elements calculation) are provided in the common block 
/cms_ef f _momdec/, as given in Table ^ They may be useful for example to impose some 
additional cuts. However, care must be taken, since such cuts might be unphysical. 

In the case of a MC run with variable-weight events, the user is provided with a main 
weight through 



CALL KW_GetWtMain(WtMain) 



! get main Monte Carlo weight 



Of course, for constant- weight runs WtMain = 1. For special purposes the user may be 
also interested in auxiliary weights, which are provided with the help of 



31 



CALL KW_GetWtAll(WtMain,WtCrud,WtSet) ! get all Monte Carlo weights 
where REAL*8 WtSet(lOO) is an array of weights described in Table ITU|. Alternatively 



these weights are available from the common block /wtset/, see Table |T0|. The aux- 
iliary weight should be defined as WtCrud*WtSet (i), and the corresponding integrated 
cross section simply obtained by multiplying its average by the crude normalization cross 
section, which is provided through 

CALL KW_GetXSecNR(XSecNR,XErrNR) ! get normalization x-section 

The complete description of post-generation output parameters from KW_Finalize is col- 
lected in Table |TT|. 

Finally, the description of input-output parameters of the semi-analytical routine 
korwan is given in Table |l^. It should be recalled here, that the semi-analytical functions 
slwan and sls2wan are initialized by the korwan routine. To this end, for pure initializa- 



tion, a special setting of its KeyMod argument, KeyMod > 10000, is provided. Tables [12 
and |l^ describe the input /output of the mavrg and mloss routines. 

4.3 Printouts of the Program 

In this section we describe a printout of the demonstration program KWdemo in the all-four- 
fermion mode, shown in Appendix C The printout starts with the detailed specification 
of the actually used input parameters. Also, logos of all the activated libraries (TAUOLA, 
PHOTOS, JETSET) are printed here. Next, the printout of one full event (in the standard 
PDG convention) is shown. The final reports of KoralW are collected in four windows: 
A, B, C and D. 

Window A provides a technical internal report of the karlud routine, i.e. information 
on the crude distribution and the Born matrix element. The line aO gives the total number 
of generated weighted events. The line al shows the number of events with the negative 
weight wtcrud. This entry should be equal to zero. The line a2 provides the value of the 
master crude distribution. The line a3 shows the average of the weight wtcrud, and the 
line a4 - the corresponding "cross section" (no matrix element is included here - only the 
crude distribution formula). The lines a5-a8 repeat the information of the lines aO-al and 
aS-aA, but for the weight with the Born matrix element, i.e. wtcrud*wtset(l). The lines 
a9-all provide information on how many of the events have the weight wtcrud*wtset(l) 
over the wtmax, the maximal weight for rejection. Their contribution to the cross section 
in absolute units (picobarns) and relative to the Born cross section is listed in alO and 
all, respectively. Note that there are no /3-functions included in any of the printouts in 
Window A. 

Window B is devoted to the technical information on the QED radiative corrections 
in different orders in a within the YFS framework. This information is quite important, 
because it shows how big the contributions of subsequent orders of the perturbative series 
are, and thus allows us to estimate the missing higher-order effects. The entries 63- 
66 provide the values of the total cross section (in picobarns) in the orders 0{a^)exp- 
0{a^)exp, while the entries 617-619 give the differences of these values in subsequent orders. 
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The entries 67-616 contain the information on the contributions from the individual YFS 
residuals /3i, also in different orders in a. The differences between various residuals and 
various orders for a given residual are provided in the entries 620-625. 

Window C is the most important from the user's point of view. It provides the total 
number of generated events (cl), the number of accepted events (c2), the best-order total 
cross section (in picobarns) with the absolute error (c3), as well as the relative error (c4), 
for the generated statistics sample. Then the information on the negative-weight events is 
given as a number of such events (c5) and their relative contribution to the cross section 
(c6), followed by similar information on the overweighted events (c7 and c8, respectively). 
All these four numbers should be as small as possible. 

Window D contains some supplementary information and is printed out only if the 
4-fermion matrix element is included in the calculation. Part I provides the values of the 
average weight idV) and the total cross section ((i2) for the WW process (if such exist for a 
given final state). Part II contains some additional information on the 4-fermion process, 
as: the average Born level weight ((i3), the average total weight ((i4), the total cross section 
((i5) - this should be the same as in Window C (entry c3), the relative contribution to 
the cross section of the non-H^l^ diagrams, i.e. the size of the background to the W^-pair 
production ((i6, calculated from entries d1 and rfS, and dl ^ calculated from a dedicated 
weight, monitored during the event generation). 

Finally, the report on different channels is printed out. It includes the codes of the 
decay W and/or Z channels (first column), the flavours of the final-state fermions in the 
"human-readable" format (second column), the values of the total cross section and its 
absolute error (in pb) for individuals channels (third column) , the ratio of the maximum 
weight over the average weight (fourth column) and over the average non-zero weight 
(fifth column), the fraction of events generated for a given channel with respect to the 
total number of events (sixth column) and a similar fraction for non- zero-weight events 
(last column). The last line of this report provides the value of the total cross section and 
its error (in pb) summed over all open channels. 

This completes the description of the output of KoralW. The remaining entries shown 
in the demo output are produced by the demo main program. 

4.4 Random Number Generators 

The KoralW code uses exclusively one of three random number generators: RANMAR, 
ECURAN or CARRAN. These single precision generators are called by one double precision 
interface routine VARRAN. The choice between generators is done with the help of the key 
KeyRnd (1 =RANMAR, 2 =ECURAN, 3 =CARRAN). In order to avoid possible interference with 
libraries, RANMAR of KoralW is renamed to MARRAN. JETSET, PHOTOS and TAUOLA 
have their own independent random number generators. 
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Appendix A: 

Program Parameters and Their Settings 



Parameter 



Position and meaning 



cmsene 

gmu 

alfwin 

amaz 

gammz 

amaw 

gammw 

vvmin 
vvmax 
wtmax 

amh 

agh 

alpha_s 

arbitr 

arbitrl 

themin 

arbitr2 

wtmax_cc03 

PReco 



xpar(l) (=180) : ^/i, centre-of-mass (CMS) energy [GeV] 

xpar(2) (=1.16639d-5) : G^, Fermi constant in [GeV] 

xpar(3) (=128.07d0): 1/aw inverse QED coupling constant at Mw scale 

xpar (4) (=91 . 1888) : M^, mass of Z boson [GeV] 

xpar(5) (=2.4974): Tz, width of Z boson in [GeV] 

xpar (6) (=80.230): Mw, mass of W boson [GeV] 

xpar (7) (=-2.03): Pv^, width of W boson [GeV], for (gammw < 0) Tw is 

recalculated from G^, Myy and as 

xpar(8) (=ld-6): Minimum t;-variable (dimensionless), infra-red cut-off 

xpar (9) (=0.99): Maximum value of f -variable 

xpar (10) (=-1): Maximum weight for rejection, for wtmax < redefined 

inside program 

xpar (11) (=1000): Higgs mass [GeV] 

xpar(12) (=1): Higgs width [GeV] 

xpar(13) (=0.12): QCD coupling constant 

xpar (14) (=600): minimum visible pT [GeV^] 

xpar (15) (=8): invariant-mass cut for e^e~ ff [GeV^] 

xpar (16) (=ld-6): minimum 9 [rad] relative to beam (0=no cut) 

xpar (17) (=300): maximum pT of photons in e^e^ff channels [GeV^] (if 

arbitr2 < then no cut) 

xpar (18) (=-1) : maximum CC03 weight for rejection in KeySmp=0, for wtmax 

< redefined inside program 

xpar (19) (=0): Colour Reconnection Probability 



Table 1: List of input parameters of the KoralW generator in xpar vector. Default values 
in brackets. 
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Parameter 


Position and meaning 




xpar (21-57): Values of TGCs: set 1, most general set - complex numbers, 




default values are wild random, not shown 


gl(l) 


=DCMPLX(xpar(21) ,xpar(31)) = gf, for WWZ vertex 


kap(l) 


=DCMPLX(xpar(22) , xpar (32)) = Kz 


for WWZ vertex 


lamd) 


=DCMPLX(xpar(23) , xpar (33)) = A^ 


for WWZ vertex 


g4(l) 


=DCMPLX(xpar(24) , xpar (34)) = g^ 


for WWZ vertex 


g5(l) 


=DCMPLX(xpar(25) ,xpar(35)) = ^if 


for WWZ vertex 


kapt(l) 


=DCMPLX(xpar(26) , xpar (36)) = k^ 


for WWZ vertex 


lamtd) 


=DCMPLX(xpar(27) , xpar (37)) = A^ 


for WWZ vertex 


gl(2) 


=DCMPLX(xpar(41),xpar(51)) = gf 


for WW J vertex 


kap(2) 


=DCMPLX(xpar(42) , xpar (52)) = Kg 


for WWj vertex 


lam(2) 


=DCMPLX(xpar(43) , xpar (53)) = A^ 


for WW J vertex 


g4(2) 


=DCMPLX(xpar(44),xpar(54)) = gl 


for WWj vertex 


g5(2) 


=DCMPLX(xpar(45),xpar(55)) = 5(f 


for WW J vertex 


kapt(2) 


=DCMPLX(xpar(46) , xpar (56)) = kg 


for WWj vertex 


laint(2) 


=DCMPLX(xpar(47) , xpar (57)) = A^ 


for WW'j vertex 




xpar(61-65): Values of TGCs: set 2, see CERN 96-01, Vol. 1, p. 525 


delta_Z 


=xpar(61) = 6z 


x_gainina 


=xpar(62) = x^ 


x_Z 


=xpar(63) = xz 


y_gainma 


=xpar(64) = y^ 


y-z 


=xpar(65) = yz 




xpar (71-73): Values of TGCs: set 3, see CERN 96-01, Vol. 1, p. 525 


alpha_Wphi 


=xpar(71) = aw(i> 


alpha Rphi 


=xpar(72) = a^^ 


alpha_W 


=xpar(73) = aw 



Table 2: List of input parameters of the KoralW generator in xpar vector (cont.). Default 
values in brackets. 
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Variable 



Position and meaning 



KeylSR 



KeyFSR 
KeyNLL 



KeyCul 
KeyBra 



KeyMas 

KeyZet 

KeySpn 
KeyRed 



KeyWu 



xpar(lOll) (=1) 
=0 Initial State Radiation is OFF 
=1 Initial State Radiation is ON 

xpar(1012) (=0) Final State Radiation switch, INACTIVE 
xpar(1013) (=1) 

=0 sets to zero the Next-to Leading a/vr terms in the YFS form factor, useful 
for comparisons 

=1 the q/tt terms are kept in the YFS form factor 
xpar(1014) (=1) 
=0 Coulomb correction is OFF 
=1 Coulomb correction is ON 
xpar(1021) (=1) 

sets W branching ratios, used for normalization of CC03 matrix element only 
=0 Born values (no mixing) with Naive QCD (if as = 0: ud, cs = 1/3, 
ezv, fiu, Tu = 1/9, others = 0) 

=1 with CKM mixing and Naive QCD, defined in KW_Initialize 
=2 with CKM mixing and Naive QCD, calculated in IB A from the CKM matrix 
(PDG '98) 
xpar(1022) (=1) 

=0 Massless kinematics; r decay, radiative corrections in decay and hadroniza- 
tion must be switched off 

=1 Massive kinematics, masses of fermions are read from xpar 
xpar(1023) (=0) 

=0 Z width in Z propagator: {s/Mz)Tz 
=1 Z width in Z propagator: MzTz 
=2 no Z width in Z propagator 
xpar(1024) (=1) 

=0 spin effects are switched OFF in W decays, for tests only 
=1 spin effects are switched ON in W decays 
xpar(1025) (=0) 

=0 "sophisticated" reduction of massive to massless four-vectors for CC03 
internal Matrix Element 

=1 "brute-force" reduction of massive to massless four-vectors for CC03 inter- 
nal Matrix Element, four-momentum NON-conserving 

=2 no reduction at all of massive to massless four-vectors for CC03 internal 
Matrix Element 
xpar(1026) (=0) 

=0 W width in W propagator: [s/MyY)Tyi/ 
=1 W width in W propagator: MyyTw 
=2 no W width in W propagator 



Table 3: List of input parameters of the KoralW generator in xpar vector (cont.). Default 
values in brackets. 
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Parameter 



Position and meaning 



KeyWgt 



KeyRnd 



KeySmp 



KeyMix 

Key4f 
KeyAcc 



KeyZon 
KeyWon 



xpar(1031) (=0) 

=0 constant weight wtmod=l, for apparatus Monte Carlo 

=1 variable weight wtmod events 

=2 for special purposes: wtmod=l for internal matrix element, AND varying 

weight for external matrix element 

xpar(1032) (=1) 

=1 RANMAR random number generator 

=2 ECURAN random number generator 

=3 CARRAN random number generator 

xpar(1033) (=2) 

=0 presampler set as in KoralW v. 1.02-1.2x, i.e. CC03 oriented 

=1 first presampler for all 4-fermion final states 

=2 second presampler for all 4-fermion final states 

=3 50/50 mixed (1+2) presampler for all 4-fermion final states 

xpar(1041) (=0) 

=0 "renormalization scheme" of LEP2 Workshop (recommended) 

=1 "renormalization scheme" based on G^ 

=2 "bare Born" (for tests) 

xpar(1042) (=1) 

=0 External Matrix Element OFF 

=1 External Matrix Element ON 

xpar(1043) (=0) 

=0 anomalous WWV couplings in internal CC03 matrix element OFF 

>0 anomalous WWV couplings in internal CC03 matrix element ON 

=1 the most general (complex number) TGCs in the notation of K. Hagiwara 

et al., Nucl. Phys. B282 (1987) 253 

=2 parametrization of CERN 96-01, Vol. 1, p. 525: 6z,x.y,xz,y'y,yz 

=3 parametrization of CERN 96-01, Vol. 1, p. 525: aw(f>,ctB(f>,Oiw 

xpar(1044) (=1) 

KeyZon=0 ZZ type final states OFF 

KeyZon=l ZZ type final states ON 

xpar(1045) (=1) 

KeyWon=0 WW type final states OFF 

KeyWon=l WW type final states ON 



Table 4: List of input parameters of the KoralW generator in xpar vector (cont.). Default 
values in brackets. 
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Parameter 



Position and meaning 



KeyDwm 



KeyDwp 



Nout 
Jakl 
Jak2 



itdkrc 



ifphot 



ifhadM 
ifhadP 



Umask 



xpar(1055) (=0) 

Sets decay channel of W" or "Zi" resonance, depending on KeyWon/KeyZon 

=0 inclusive, otherwise exclusive modes for W are 



= 12 3 


4 


5 


6 


7 


ud cd us 


cs 


ub 


cb 


ez^ 


and for Z are 










= 12 3 


4 


5 


6 


7 


dd uu ss 


cc 


bb 


ee 


/x^ 



flV TV 



TT 



= 9 



10 



11 

VtUt 



xpar(1056) (=0) 

Sets decay channel of W~^ or "Z2" resonance, depending on KeyWon/KeyZon, 
similar assignments as for W~ or "Zi" 

Note: In the inclusive assignment KeyDwm=0, KeyDwp=0 the auxiliary user 
mask Umask is used by the program. It allows the user to config- 
ure arbitrary menu of KeyDwm, KeyDwp. The Umask is passed through 
xpar ( 1101) -xpar (1302) parameters. Final states that may come from ei- 
ther WW or ZZ are classified as VFPF-type except for ussu, ubbu, cddc and 
cbbc, which are classified as ZZ-type. See Appendix B for more details. 
xpar(1057) (=-1) 

Output unit number for the generator (if < then Nout=16) 
xpar(1071) (=0) 

Input for TAUOLA, defines decay mode of r^ in W^ decay 
xpar(1072) (=0) 

Input for TAUOLA, defines decay mode of t~ in W~ decay 
Jakl, Jak2 = -1 TAUOLA is switched OFF 
Jakl , Jak2 = requests all r decay channels to be simulated 
Jakl , Jak2 > single specific r^ decay channel, see TAUOLA manual 
xpar(1073) (=1) 

Input for TAUOLA, radiative corrections in leptonic r decays switch 
itdkrc=l corrections are ON 
itdkrc=0 corrections are OFF 
xpar(1074) (=1) : PHOTOS, activation switch 
=1 PHOTOS is ON 
=0 PHOTOS is OFF 

xpar(1075) (=1) : W~ hadronization activation switch (JETSET) 
xpar(1075) (=1) : W^ hadronization activation switch (JETSET) 
ifhadm,ifhadp=l hadronization is ON 
ifhadm,ifhadp=0 hadronization is OFF 
In the present version if hadm and ifhadp have to be equal! 
xpar(1101-1302) 
user supplied Umask to activate a specific menu of the final states. 



Table 5: List of input parameters of the KoralW generator in xpar vector (cont.). Default 
values in brackets. 
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Parameter 



Meaning 



amel 

alf inv 

gpicob 
br(l:20) 

amafin(20) 



vckm(l:3,l:3) 



xpar(lOO) (=0.51099906d-3) Electron mass used by the bremsstrahlung 

generator, reset to amaf in(ll) 

xpar(lOl) (=137.0359895) l/a(0) Inverse QED coupling constant at 

Thomson scale Q^ =0 

xpar(102) (=389.37966d6) Conversion factor from GeV"^ to pb 

xpar(131-139) W branching ratios; numbering of entries is: 1 = ud, 2 = 

cd, 3 = us, 4 = cs, 5 = u6, 6 = c5, 7 = e, 8 = /U, 9 = r 

xpar(500 +10*KF +6) Masses of the W decay products; used entries KF 

(in PDG notation) are: 1 = d, 2 = u, 3 = s, 4 = c, 5 = 6, 6 = t, 

11 = e, 12 = zvg) 13 = /i, 14 = f^, 15 = r, 16 = Vt\ masses of r and v^- have to 

be independently set to the same numerical values in intialization of TAUOLA 

xpar( 111-119) CKM matrix elements 



Table 6: List of input parameters of the KoralW generator in xpar vector (cont.). Default 
values in brackets. 



Parameter 



Meaning 



qeffl(4) 

qeff2(4) 

sphum(4) 

sphot(100,4) 

nphot 



Effective parameter for matrix element (e~ effective beam, CMS) 
Effective parameter for matrix element (e"*" effective beam, CMS) 
Sum of four-momenta of ISR photons, CMS 
Four-momenta of ISR photons, CMS 
Multiplicity of ISR photons 



Table 7: List of four-momenta in the internal common block /momset/ of the KoralW 
generator. 



Parameter 


Meaning 


ql(4) 


Four-momentum of the W^ /Zi resonance, = pi +P2, in CMS 


q2(4) 


Four-momentum of the W^ /Z2 resonance, = Ps + P4, in CMS 


pi (4) 


Four-momentum of the fermion from W~ /Zi decay, CMS 


p2(4) 


Four-momentum of the antifermion from W~ /Zi decay, CMS 


p3(4) 


Four-momentum of the fermion from W^ /Z2 decay, CMS 


p4(4) 


Four-momentum of the antifermion from W~^ /Z2 decay, CMS 



Table 8: List of four-momenta in the internal common block /momdec/ of the KoralW 
generator. 
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Parameter 


Meaning 


ef f beaml (4) 


Four-momentum of the e~ beam, along z-axis, in effective CMS 


effbeam2(4) 


Four-momentum of the e"*" beam, in effective CMS 


effpl(4) 


Four-momentum of the fermion from W~ /Zi decay, effective CMS 


effp2(4) 


Four-momentum of the antifermion from W^ /Zi decay, effective CMS 


effp3(4) 


Four-momentum of the fermion from W^ /Z2 decay, effective CMS 


effp4(4) 


Four-momentum of the antifermion from W^ /Z2 decay, effective CMS 



Table 9: List of four-momenta in the internal common block /cms_eff jnomdec/ of the 
KoralW generator. 



Parameter 



wtcrud 
wtmod 



wtset(l-lOO) 



wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 
wtset 



(1) 

(2) 

(3) 

(4) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(40) 

(41-49) 



Meaning 



Crude weight, necessary to build the total weight out of wtset 

Best weight 

For KeyISR=0 wtmod=wtcrud*wtset(l). 

For KeyISR=l wtmod=wtcrud*wtset (4) 

Born matrix element with various 13 contributions - to get total 

weight must be multiplied by wtcrud 

Zero-order complete (/3o) 

First-order complete (/3o + Pi) 

Second-order complete (^^o + f3i + (32) 

Third-order complete {f3o + I3i + (32 + (3^) 



0{a^ 
0(a^ 



contribution 

0-real, 1-virtual photon contribution 

1-real, 0-virtual photon contribution 

0-real, 2-virtual photon contribution 

1-real, 1-virtual photon contribution 

2-real, 0-virtual photon contribution 

0-real, 3-virtual photon contribution 

1-real, 2-virtual photon contribution 

2-real, 1-virtual photon contribution 

3-real, 0-virtual photon contribution 
wtniod4f , external matrix element 
wt4f(l-9) additional weights from ext. matrix element, if provided 



Table 10: List of output weights in the common block /wgtall/ of the KoralW generator. 
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Parameter 


Meaning 


XSecMC 


Principal best total Monte Carlo cross section [pb] 




For KeyISR=0 Born cross section 




For KeyISR=l third-order exponentiated cross section 


XErrMC 


Its absolute error [pb] 


XSecNR 


Normalization cross section: 




For KeyWgt=0 principal cross section XSecMC [pb], 




For KeyWgt=l crude cross section XCrude [pb] 


XErrNR 


Its absolute error [pb] 


NevMC 


Total number of generated Monte Carlo events 



Table 11: List of output parameters of KoralW generator transferred through parameters 
of the getter type subprograms KW_GetXSecMC (XSecMC , XErrMC) , KW_GetNevMC (NevMC) and 
KW_GetXSecNR(XSecNR, XErrNR) . 



Parameter 


Meaning 




INPUT 


svar 

keypho 

kaccbre 


s, CMS energy squared [GeV^] 
As KeyMod of korwan 
As KeyPre of korwan 




OUTPUT 


dmavrg 
dmerr 


Mass average in GeV 
Absolute error in GeV 



Table 12: List of arguments of the mavrg routine. 



Parameter 


Meaning 




INPUT 


svar 

keypho 

eeps 


s, CMS energy squared [GeV^] 

As KeyMod of korwan 

Absolute precision of integration 




OUTPUT 


vvloss 
vverr 


Average mass loss in GeV 
Absolute error in GeV 



Table 13: List of arguments of the mloss routine. 
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Parameter 


Meaning 




INPUT 


svar 
vvmin 
vvmax 
keymod 

keypre 


■s, CMS energy squared [GeV'^] 

Vmim minimal v variable, in most cases should be set to 

Vmax, maximal v variable 

Defines type of structure functions used for ISR: 

= No ISR, Born, 

=300 Zero Order, YFS style, 

=301 First Order, YFS style, 

=302 Second Order, YFS style, 

=303 Third Order, YFS style, 

=502 Second Order, Gribov-Kuraev-Fadin style, 

=310 First-Order YFS BetaO only, 

=311 First-Order YFS Betal only, 

=320 Second-Order YFS BetaO only, 

=321 Second-Order YFS Betal only, 

=322 Second-Order YFS Beta2 only, 

<0 as (-keymod) but multiplied by v differential distribution da/dlogv 

Defines precision level of the computation 

For KeyMod=0 No ISR, Born, in (ez^e) channel 

=1 absolute error 1 x 10^ [pb] 

=2 absolute error 1 x 10^^ [pb] 

=3 absolute error 1 x 10^^ [pb] 

For KeyMod>0 ISR, in (ez^e) channel 

=1 absolute error 3 x 10^ [pb] 

=2 absolute error 1 x 10"^ [pb] 

=3 absolute error 1 x 10^^ [pb] 

=4 absolute error 1 x 10"'^ [pb] 




OUTPUT 


xsect 
errabs 


Cross section pb 
Absolute error [pb] 



Table 14: List of arguments of the korwan routine. 
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Appendix B: 

Umask matrix for Final States— Default values 

In this appendix we present the template mask for specifying final states in the fully inclu- 
sive A-fermion mode, i.e. with KeyWon=KeyZon=l, KeyDwm=KeyDwp=0. The mask below 
opens all channels allowed by the program. To exclude some of them, appropriate zeros 
should be put into the Umask matrix. There is a total of 202 entries in it, corresponding to 
an easy-to-handle classification: 9 x9 oiW xW decay modes and 11 xll oi Z x Z decay 
modes. This way some decay channels are doubly counted as some final states can be both 
WW- and ZZ-tjpe. To avoid such double-counting in the event generation, appropriate 
entries in the mask are set to 0. Regardless of the input values used, the program always 
checks for these MIX-type channels and always assigns them either to WW or to ZZ. 
Therefore the user can safely use the mask with all entries set to 1. The actual mask used 
by the program after this internal verification is printed to the output file. Note also that 
we adopted a convention that four of the MIX-type final states {ussu, ubbu, cddc, cbbc) are 
coded as ZZ (xpar(1205), xpar(1215), xpar(1227) and xpar(1229)) whereas other 
MIX-types as WW (cf. Subsection p.2.2|) . Finally, it should be kept in mind that this is 
one mask for all channels and not two separate masks for WW and ZZ ones. 



xpar (1101-1302)= Umask: 



Wm= 1 


:ud 


2:cd 


3: us 


4:cs 


5 rub 


6:cb 


7: el 


8:mu 


9:ta 






/Wp= 


xpar(1101-1109) : 




1 


1 


1 




1 


1 










l:ud 


xpar(1110-1118) : 







1 


1 




1 


1 










2:cd 


. . . 




1 





1 




1 


1 










3:us 






1 


1 


1 




1 


1 










4:cs 






1 


1 


1 




1 


1 










5:ub 






1 


1 


1 







1 










6:cb 






1 


1 


1 




1 


1 










7:el 






1 


1 


1 




1 


1 










8:mu 


xpar(1173-1181): 




1 


1 


1 




1 


1 










9:ta 


Zl= 1 


:d 


2:u 


3:s 


4:c 


5:b 


6: el 


7:mu 


8:ta 


9:ve 


lOvm 


llvt 


/Z2= 


xpar(1182-1192) : 


































l:d 


xpar(1193-1203) : 




1 





























2:u 


. . . 




1 


1 


























3:s 






1 





1 























4:c 






1 




1 






















5:b 






1 




1 




1 

















6: el 






1 




1 




1 


1 














7:mu 






1 




1 




1 


1 


1 











8:ta 






1 




1 







1 


1 


1 








9:ve 






1 




1 




1 





1 


1 


1 





lOvm 


xpar(1292-1302): 




1 




1 




1 


1 





1 


1 


1 


llvt 
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### ### 


* 


## ### ### 


* 


## ### ### 


* 


## ### ### 


* 


## ### # ### 


* 


## ### ### ### 


* 


####### #### #### 


* 


## ## 


* 



Appendix C: Output of the Demo Program 

* ### ### #### ###### ## 

* ### ### ## ## ## ## #### 

* ###### ## ## ## ## ## ## 

* ###### ## ## ##### ## ## 

* ### ### ## ## ## ## ######## 

* ### ### #### ## ## ## ## 

* ### ### version 1.42.3 

*^^^^^^^^^^^^^^^^^^^^^ ivl ^ "Y^ /^r^ n Q Q Q ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

Last modification: 3.17.1999 

* Written by: * 

* S. Jadach (Stanislaw.Jadach@cern.ch) * 

* W. Placzek (Wieslaw.Placzek@cern.ch) * 

* M. Skrzypek (Maciej .Skrzypek@cern. ch) * 

* B.F.L. Ward (bflw@slac.stanford.edu) * 

* Z. Was (Zbigniew.Was@cern. ch) * 

* Papers : * 

* M. Skrzypek, S. Jadach, W. Placzek, Z. Was * 

* CERN-TH/95-205, Jul 1995, CPC 94 (1996) 216 * 

* M. Skrzypek, S. Jadach, M. Martinez, W. Placzek, Z. Was * 

* CERN-TH/95-246, Sep 1995, Phys . Lett. B372 (1996) 289 * 

* S. Jadach, W. Placzek, M. Skrzypek, B.F.L. Ward, Z. Was * 

* CERN-TH/98-242, UTHEP-98-0702, Jul 1998, submitted to CPC * 

* M. Skrzypek, S. Jadach, W. Placzek, B.F.L. Ward, Z. Was * 

* CERN-TH/99-06, UTHEP-98-1001 , Jan 1999, proc. of RADC0R98 * 

* Related papers: * 

* T. Ishikawa, Y. Kurihara, M. Skrzypek, Z. Was * 

* CERN-TH/97-11, Jan 1997, Eur. Phys. J. C4 (1998) 75 * 

* S. Jadach, K. Zalewski * 

* CERN-TH/97-29, Jan 1997, Acta Phys. Pol. B28 (1997) 1363 * 

* WWW : * 

* \protect\vrule widthOpt\protect\href {http://hpjmiady. if j .edu.pl/}{http://hpjm: 

* Acknowledgements: * 

* We acknowledge warmly very useful help of : * 

* M. Martinez in testing versions 1.01 and 1.02, * 

* M. Gruenewald and A. Valassi in testing version 1.21 * 

* S. Jezequel in testing versions 1.31-1.33 * 

* M. Witek in testing version 1.41 * 

* M. Verzocchi in testing version 1.42 * 

* KORALW input parameters used * 

* 190.00000000 CMS energy total CMSENE I.O * 

* 1101 QED super-switch KeyRad IQl * 

* 1 Init. state Rad. KeylSR IQ2 * 

* Final state Rad. KeyFSR IQ3 * 

* 1 Next. To Leading KeyNLL IQ4 * 

* 1 Coulomb corr. KeyCul IQ5 * 

* 1012 Physics super-switc KeyPhy IPl * 

* FS mass reduction KeyRed IP2 * 

* 1 Spin in W decays KeySpn IP3 * 

* Z propag. KeyZet IP4 * 

* 1 Mass kinematics. KeyMas IP5 * 

* 2 Branching Rat. KeyBra IP6 

* W propag. KeyWu IP7 * 
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* 



* 
* 
* 
* 
* 
* 

* 
* 
* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 
* 
* 

* 
* 
* 
* 

* 
* 
* 
* 
* 

* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 



211 Technical super-swi 

2 presampler type 

1 rand Numb type 

1 weighting switch 

11010 Miscelaneous 

sinW2 input type 

1 4 fermion matr el 

Anomalous couplings 

1 WW type final state 
1 ZZ type final state 

W-/Z decay mode 

W+/Z decay mode 

1.16639000 G_mu * ld5 

128.07000000 inv alpha_w 

91.18880000 Z mass [GeV] 

2.49740000 Z width [GeV] 

80.23000000 W mass [GeV] 

2.08545732 W width [GeV] 

.00000100 dummy infrared cut 

.99000000 v_max ( =1 ) 

2.00000000 max wt for rejectn. 

4.00000000 max wt for CC03 rej 

.12000000 alpha_s: QCD coupl. 

.00000000 Color Re-Con. Prob. 

.23103091 sin(theta_W)**2 

Z width in Z propagator: s/M_Z *GAMM_Z 

CKM matrix elements: 

.97525000 V_ud 

skipped .... 

.99925000 V_tb 
Unitarity check of the CKM matrix: 

1.000 .000 .003 

VV+ = .000 1.000 .000 

.003 .000 1.000 



* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 



.32097393 

skipped . . . 

.10835195 

.01000000 

.00500000 

skipped . . . 

1.77710000 
.00100000 



Branching ratios of W decays: 



fermion masses: 



ud 

tau 

d 

u 

tau 
vtau 



KeyTek 
KeySmp 
KeyRnd 
KeyWgt 

KeyMis 

KeyMix 

Key4f 

KeyAcc 

KeyWon 

KeyZon 

KEYDWM 
KEYDWP 

GMU 

ALFWIN 

AMAZ 

GAMMZ 

AMAW 

GAMMW 

VVMIN 

VVMAX 

WTMAX 

WTMAX 

ALPHAS 

PReco 

SINW2 



VCKM(1,1) 
VCKM(3,3) 



* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 



Predefined cuts on final state fermions 

600.00000000 min. vis p_t**2 

8.00000000 add. cut for e+e-ch+ 

, lOOOOOOOE-05 min. theta with beam 

300.00000000 max. p_t**2 phot eex 

DECAY LIBRARIES 

TAUDLA for W+ 

TAUDLA for W- 

1 TAUDLA Drd( alpha) 
1 PHDTDS 
1 JETSET for W- 
1 JETSET for W+ 



BR(1) 
BR(9) 

AMAFIN(l) 
AMAFIN(2) 

AMAFIN(15) 
AMAFIN(16) 

GeV~2 
GeV~2 

rad 
GeV-2 

JAKl 
JAK2 
ITDKRC 
IFPHDT 
IFHADM 
IFHADP 



ITl * 
IT2 * 
IT3 * 
IT4 * 

IMl * 
IM2 * 
IM3 * 
IM4 
IM5 
IM6 



IDl 
ID2 



* 
* 
* 
* 

* 
* 
* 
* 

* 
* 
* 
* 

* 
* 

* 

* 
* 
* 

* 
* 

* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 

IBl * 
IB9 * 

IMl * 
IM2 * 



1.9 
110 
111 
112 

113 



IVl 



IV9 



IMIO 
IMll 

X2 
X3 
X6 
X3 

ILl 
IL2 
IL3 
IL4 
IL5 
IL6 



umask_init=>umask : 
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 
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1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 



.0 
1.0 
















_ 
1.0 
1.0 



1.0 

.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.0 

.0 

1.0 

.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 



1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
.0 
.0 
.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 



1.0 
1.0 
1.0 



1.0 

1.0 

.0 

.0 

.0 

.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 



1.0 

.0 

1.0 

1.0 

1.0 

.0 

.0 

.0 

.0 

.0 

1.0 

1.0 

1.0 

.0 

1.0 



1.0 
1.0 
1.0 



1.0 

1.0 

.0 

.0 

.0 

.0 

.0 

.0 

1.0 

1.0 

1.0 



1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

1.0 

1.0 



.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
1.0 



.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 



1.0 1.0 1.0 1.0 1.0 1.0 1.0 .0 1.0 1.0 1.0 



* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

* 
* 
* 
* 
* 



* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

JAKl = DECAY MODE FERMIONl (TAU+) * 

JAK2 = DECAY MODE FERMI0N2 CTAU-) * 



*****TAUOLA LIBRARY: VERSION 2.6 ****** 
***********August 1995*************** 
**AUTHORS: S.JADACH, Z. WAS************* 
**R. DECKER, M. JEZABEK, J .H.KUEHN***** 
**AVAILABLE FROM: WASM AT CERNVM ****** 
***** PUBLISHED IN COMP . PHYS . COMM.*** 
*******CERN-TH-5856 SEPTEMBER 1990***** 
*******CERN-TH-6195 SEPTEMBER 1991***** 
*******CERN-TH-6793 NOVEMBER 1992***** 
**5 or more pi dec: precision limited 
******DEXAY ROUTINE: INITIALIZATION**** 
JAKl = DECAY MODE FERMIONl (TAU+) 
JAK2 = DECAY MODE FERMI0N2 (TAU-) 



Window H used only by Grace 2.0 * 

Higgs boson parameters * 

1000.00000000 xpar(ll)= higgs mass amh HI * 

1.00000000 xpar(12)= higgs widt agh H2 * 

* 



* Window X_ZZ * 

* mm_brancher_ZZ report * 

* mm_brancher_ZZ is on * 

* .00000000 Drob. for branch NR: 1 XI * 



.00000000 
skipped . . . 

.02232143 
skipped . . . 

.00000000 



prob. for branch NR 
prob. for branch NR 
prob. for branch NR 



1 
8 
65 



XI 
XI * 
XI * 



* 



Event list ing ( st ajidard) 
I particle/ jet K(l,l) K(l,2) K(l,3) KCl,4) 



Czo) 




11 


Czo) 




11 


e- 




1 


e+ 




1 


mu- 




1 


mu+ 




1 




sum 


charge 



sum momentum and inv. mass: 





00000 




00000 


95.00000 


95.00000 




00061 




00000 




00000 


-95 . 00000 


95.00000 




OOOBl 


41 


18296 


13 


10538 


23.92689 


129.21960 


119 


40434 


41 


18296 


-13 


10B38 


-23 . 92689 


60 . 78049 


35 


41166 




01444 




00373 


-46.54470 


46.54470 




00061 


41 


19740 


13 


10166 


70.47158 


82.67480 




OOOBl 


-3 


64304 


1 


49896 


4.29885 


5.83179 




10666 


44 


82600 


-14 


60434 


-28 . 22573 


54 . 94870 




10666 




00000 




00000 


.00000 


190.00000 


190 


00000 



* Window X_WW * 

* mm_brancher_WW report * 

* mm_brancher_WW is on * 

* .46296296 prob. for branch NR: 1 XI * 
skipped .... 
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* .00000000 prob. for branch NR: 65 XI * 



KORALW final report 

Window A 
WEIGHTED evts. 

ecru matrix element means : 

a) Born matrix element for CC03 processes 

b) technical crude m.e. for nc processes or 
for keysmp .NE. 



xsect with no matrix element 



10000 

35058.681 
.55820244 
19569.841 



xsect with ecru matrix el 
10000 



.40047631E-03 
14.040171 



total no of events 
wtcrud < evts 
sigma_crude 
+- .64458 <wtcrud>, rel err 
+- 12614. phsp. vol, no beta-0 

only, no betas 
total no of events 
wtcrud* wtborn <0 evt 
+- .27055 <wtcrud*wtborn>, rel 
+- 3.7986 sigma (born m.el.) 



* 
* 
* 
* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

* KDRALW final reoort * 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 
* 
* 
* 
* 
* 
* 
* 



.OOOOOOOOE+00 
.OOOOOOOOE+00 



xsect over wtmax_cc03 
ecru matrix el. only, no betas 

evts: wt>wtmax_ee03 
+- .OOOOOE+00 sigma: wt>wtmax_ee03 
+- .OOOOOE+00 relat sigma: wt>wtma 







* 






* 






* 






* 






* 






* 






* 






* 


nevtot 


aO 


* 


nevneg 


al 


* 


Xerude 


a2 


* 


wtkacr 


a3 


* 


xskr 


a4 


* 
* 
* 
* 


nevtot 


a5 


* 


nevneg 


a6 


* 


wtkabo 


a7 


* 


xskaO 


a8 


* 
* 
* 
* 


nevove 


a9 


* 


xskabo 


alO 


* 


xskabo 


all 


* 



KORALW final report 
Window B 
Xsec-s in [pb] 



29.68888933 

30.29512756 

30.33902996 

30.33978251 









29.68888933 

31.38896196 

. skipped . . . 

.02364120 
-.00003663 

.60623824 
.04390240 
.00075255 



+- 
+- 
+- 
+- 



+- 
+- 

+- 
+- 

+- 
+- 
+- 



4.74480757 
4.95372145 
4.96038178 
4.96052513 



4.74480757 
5.01650913 

.01578300 
-.00003099 



xsee total 
xsee total 
xsee total 
xsee total 
wt<0 events 
wt<0 events 
wt<0 events 
wt<0 events 
xsee(betaOO) 
xsee(betaOl) 

xsee(beta21) 
xsee(beta30) 



xsee_tot differences 
,29779330 xstot (alf 1-0) 
.01499107 xstot (alf 2-1) 
.00022070 xstot (alf 3-2) 

betas differences 
.27170156 xs(betaOl-OO) 



1.70007263 +■ 
skipped .... 

* -.00055056 +- -.00019179 

* -.00003663 +- -.00003099 



xs(beta21-20) 
xs(beta30) 







* 






* 


O(alfO) 


b3 


* 


O(alfl) 


b4 


* 


0(alf2) 


b5 


* 


0(alf3) 


b6 


* 


O(alfO) 




* 


O(alfl) 




* 


0(alf2) 




* 


0(alf3) 




* 


O(alfO) 


b7 


* 


O(alfl) 


b8 


* 


0(alf3) 


bl5 


* 


0(alf3) 


bl6 


* 


O(alfl) 


bl7 


* 


0(alf2) 


bl8 


* 


0(alf3) 


bl9 


* 


O(alfl) 


b20 


* 


0(alf3) 


b24 


* 


0(alf3) 


b25 


* 






* 
* 
* 
* 
* 



10000 
10000 

30.339783 



KORALW final report 
Window C 

BEST order total xsect. 

total no of events 
accepted events 

4.9605 sigma_tot [pb] 



nevtot 
NevTru 

xskabo 



* 

* 

* 

el * 

e2 * 

e3 * 
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* .16349903 

* 

* .OOOOOOOOE+00 

* 

* .OOOOOOOOE+00 



relative error errela c4 * 

events: wt<0 nevneg c5 * 

.OOOOOE+00 xsec/xtot: wt<0 xsneg c6 * 

events: wt>wtmax nevove c7 * 

.OOOOOE+00 xsec/xtot: wt>wtniax xsove c8 * 






* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 



KDRALW final report 
Window D 

Complete 4-fermion process 

I. Best ord. W-pair total xsect . 
.41443001E-03 +- .27359 <wttww> : WW weight 
14.529370 +- 3.9751 sigma_WW, best [pb] 

II. Best ord. 4-fermion total xsect. 
21819.562 +- .74956 <wtbo4f>, rel err 
.86540000E-03 +- .16350 <wttot>,rel err 
30.339783 +- 4.9605 sigma_4f, best [pb] 
.52111161 +- .15263 sigma l-Wpair/4f erm 
.52111161 +- .12996 sigma l-Wpair/4f erm 



* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

Decay Report on Different Channels 

wt_max wt_max nev_ch nev_non0 
human sigma [pb] +- abs_err 



averwt 
xskabo 

averwt 
averwt 
xskabo 
l-d2/5 
wtbgr 



dl 
d2 

d3 
d4 
d5 
d6 
d7 



wm wp 

1 1 

2 1 

3 1 

4 1 

11 11 



dq~uq uq~dq 
dq~cq uq~dq 
sq~uq uq~dq 
sq~cq uq~dq 
skipped . . . . 
nt~nt nt~nt 



.3863394E+01+- 
.0000000E+00+- 
.0000000E+00+- 
.6681917E-09+- 



.25E+01 
.OOE+00 
.OOE+00 
.67E-09 



<wt> <wt_non0> 
.67E+01 .37E+01 



tot nev_ch 
.0011 .5455 



.67E+01 
.67E+01 
.20E+01 



.37E+01 
.37E+01 
.lOE+01 



.0000 
.0000 
.0002 



.0000 
.0000 
.5000 



.0000000E+00+- .OOE+00 . lOE+01 . lOE+01 .0000 .0000 



total xsection = 30.3397825138152157 +- 4.96052512935006362 [pb] 



30.33978251 +- 4.96052513 



== demo ============ 

MC Best, XPAR, KoralW 
End demo ========== 
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